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ABSTRACT 


Two arrays of magnetometers were operated in 
southern Africa in 1971 and 1972 to supplement deep 
Schlumberger resistivity studies. The first two- 
dimensional array of 24 three-component magnetometers 
was operated in central South Africa during Sentember 
and October, 1971. Three substorm sequences have been 
analysed. Maanetoarams and maps of Fourier transform 
amplitudes and phases show a large anomaly in the vertical 
and silos horizontal variation fields over the period 
range 24-293 minutes. The vertical field has a maximum 
near Beaufort West and the horizontal field increases 
southward to a presumed maximum just south of the array. 
The principal result is the discovery of a major conductive 
Structure under the Cape Fold Belt and deep Karroo basin. 
Transfer functions between mean horizontal field components 
and the vertical component show that induction is principally 
in phase with the normal field. The effect of the Karroo 
sediments is seen in the quadrature phase induction 
anomaly. For this and other reasons the conductive body 
under the Cape Fold Belt and southernmost Karroo is 
provisionally placed in the upver mantle. A review of 
other geophysical parameters in the area indicates that 


the induction anomaly can be associated with either a 
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Lransa tron fromers00Gtorw2,600 Myreold Tithosphere® orea 
ridge on the more conductive mantle. More quantitative 
interpretation is not possible without further observations 
SOuUtNT Of ethealO7dmarray. 

Since the variation fields at the northernmost 
stations in this array contain only small anomalous 
components, these data were used for one-dimensional 
modelling. Although the data have some properties that 
indicate a deviation from a one-dimensional situation, the 
results suaqagest a slower increase in conductivity with 
depth in the range 209-800 km in the Kaapvaal shield than 
in the global model of Banks (1972). 

The second array of magnetometers was operated 
in South-West Africa, Botswana and north-western Rhodesia 
from December 1971 to February 1972. Results from three 
periods of geomagnetic disturbance are presented. Magneto- 
grams and maps of Fourier transform amplitudes and phases 
show an anomaly in the vertical and horizontal variation 
fields over the period range 21-171 minutes. A well-defined 
reversal in the vertical field is observed between QOutjo 
and Sukses on the western edge of the array. Transfer 
functions and vertical field phase maps show a zone of 
induced currents that continues the south-westward trend 
of the Luangwa-middle Zambezi Rift, an extension of the 


East African Rift system. The conductor bends south of 
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the Okavango Delta to run in line with the Walvis Ridge 
nearly due west across South-West Africa. A study of 

the induction vectors and the half-widtns of normalized 
anomalous profiles indicates a lithospheric conductor. 
The eastern branch of the conductor parallels the known 
faults in that area. Seismicity suggests that the 
Luanqwa-middle Zambezi Rift continues to the south-west 
into Botswana along the conductivity anomaly. The west- 
ward arm of the conductor is associated with seismic 
activity. The zone of induced currents coincides with 
linear positive Bouguer anomalies. Overall geophysical 
and geological evidence indicates that the high conductivity 
may be associated with a fracture zone in the lithosphere 
and may trace an extension of the African Rift system 


along oid weak zones in the continent. 
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CHAPTER I 
BACKGROUND TO THE INVESTIGATION 


1.1 Induction theory for geomagnetic induction studies 


The time varying magnetic fields used in geo- 
magnetic induction studies such as magnetotelluric sounding, 
geomagnetic deep sounding and magnetometer array studies, 
have their origin outside the earth and diffuse through 
the earth to induce currents and produce secondary fields. 

The mathematical formulation of the necessary 
induction theory begins with Maxwell's equations. In 


MKSA units these are: 


fon Bs oo City) 
gon Aeergea 20 (1.2) 
eo Hao (es) 
oa 6 (1.4) 


In an isotropic medium: 


al 

eg (1.5) 
B = ul (leon) 
D= cE Ghar 
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With B the magnetic flux density (Wb /m<, Tesla) 
H the magnetic field intensity (A/m) 
E the electric field strength (V/m) 
BD the electric displacement (C/m*) 
0 = a the resistivity (ohm.m, &-m) 
o the electrical conductivity (Siemens/m) 
J the current density (A/mo) 
e the permittivity (F/m) 


u the magnetic permeability (H/m) 


For the period range and conductivities of 
interest in magnetometer array studies the displacement 
current a is negligible compared with J (Price, 1967). 

When the curl of eq. 1.2 is taken and eq. 1.1 and the 
Wer GLONSHO fe @dS sli ose OmeanG Sliefeuis.ed, E can be eliminated 
to reduce Maxwell's equations to the induction (or diffusion) 


equation 
—> 
vis eS (153) 


Similarly, by eliminating H the diffusion equation in E 


can be written as 


es (ye (1.9) 


Tozer (1959) showed that deep in the earth yu is close to 


. (i.e. 4n > 1077 H/m), which means the diffusion equation 
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depends only on the period of the field and the resistivity. 
The atmosphere is assumed non-conducting, with p thus 


infina te So tivat 
V H = Q Cie tek) 


above the surface of the earth. Above the ground the 


magnetic field intensity Hi can therefore be expressed as 
fas WO Caan) 


where U is a scalar potential function obeying Laplace's 


equation 
VU = 0 (Gijeaiezs) 


Within the earth, where o # 0, the propagation 
of the field is controlled by another equation. Assuming 
a time dependence of exp(iwt) with w the angular frequency 
of the inducing field and neglecting displacement currents, 


SCSle lerlec. al OenulieD scaly Des COMDIned. BO, give 


iwntl inaeeen 
r—~ - (ov(=)- 9 x A) (1.13) 


This represents the basic equation for the variation 
field within the earth. 

Schmucker (1970a) separates eq. 1.13 into two 
Darts thes nirst being the “normal “solution for variations 
above a horizontally stratified medium, and the second 


being the "anomalous" variation field superimposed on the 
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normal field due to lateral variations of conductivity :in 
the substratum. The two solutions are additive since 

Cd Isls ls aebinear equation. By writing Ho Ps for the 
anomalous parts and A and Py for the normal parts where 


joc Tat aoe alg Oo oe dtd: Deen Thamind eth at 


n ai n a 
Os ae) only, the normal equation becomes 
> 
~ TwuH 
late cs (1.14) 
Pn 


The equation for the anomalous parts becomes 


vi + (vee vx fy = jou A, +o A) (1.15) 

The normal equation is generally soluble in 
cartesian and polar coordinates, as will be shown in 
Chapter 5. The anomalous equation is, however, only 
readily soluble for configurations of high symmetry or for 
two-dimensional cases (Hobbs, 1975). 

The diffusion of the primary field into the 


Earth is related to the skin depth given by the equation 


6 | ae (Giizalios) 


This defines the depth beneath the plane surface of a 
uniform conducting half-space where the amplitude of the 
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depth decreases as the resistivity and/or period decreases. 
Skin) depths of interest in array studies vary from just 
more=than 6°Kkme= Tora tield of period’ 10 min incident on 
seawater with a resistivity of .25 ohm.m, to more than 

7 ,VOORKie fOr period 24 nours andaresistivi ty -1,000-ohm.m 
(e.g. resistivity of lower crust and uppermost mantle in 


stable continental areas). 


Te2s Conductivity. distribution sin, the earth 


I22u = Globaredeomagneticesoundind: -aesitudy pot. thnemYrada all 

Time-varying magnetic fields arising outside the 
earth which acts like a conducting body, induce electric 
currents in the earth so that magnetic fields with an 
internal origin are produced. The theory of electromag- 
netic induction incorporated in the relationship between 
the parts of transient geomagnetic variations of external 
and internal origins, as obtained by analyses of world-wide 
data, provides a means of inferring the electrical conduc- 
cavity within the earth. 

The theories involved in global electromagnetic 
induction are from Lahiri and Price (1939) and Eckhardt 
(1968) and have been reviewed by Bailey (1970). In the air 
above the earth, the magnetic field can be represented as 


-VW where W is a magnetic scalar potential satisfying 
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Laplace's equation. The magnetic scalar potential can 
be represented in terms of spherical harmonics which are 
solutions of Laplace's equation. Any solution for W can 
thus be expanded as 

+n 


W(r,8,6,t) = 2 afen(t)(Z)"+en(t) (5) ""~!Jeq(cose) eM 
n=1 m=-n 


let) 


with n the harmonic degree, m the order, 6 the colatitude, 
o the longitude, a the earth's radius and en and oe are the 
amplitudes of the external and internal parts of the field 
respectively. pi (cose) is the associated Legendre function 
of degree n and order m. 

In the conductive earth use can be made of the 


magnetic vector potential such that 
B= v xA (ire) 


Lahiri and Price (1939) have shown that where the 
GONGUGTRVELY =hSiea TKUNCTIION Of cadiuUs “only, A must be of 


the form 
A =/Y x VU Brow 


with U a scalar potential. 


Taking thesctr? of "eq. 1.8 >and using the relations 


B = uA, aoe nO andere = = aA together with the assumption 
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that free charges are not present within the conductor 


é R . i 


chissinduction.equatiron which ts-alsousatisfied by U 15 
the central equation of geomagnetic deep sounding. 
insides thesconductor, so1ution of eq. 1220°for 


U yields 


Mr) p™(cose)e'™? Gla 21) 


where ea Satisfies the radial induction equation 


SM) = [n(n#1) + Gonjo(r)r@]F™ (1.22) 


Under the boundary conditions that the tangential 
components of E and H and the normal component of B are 
continuous at the surface of the sphere, the potential and 
its radial derivatives above and below the surface of the 
Sphere are fitted smoothly together. The resulting ratio 
of the induced amplitude ee toethe inducing fleld ampli tude 
en is Ss (w) the geomagnetic response of the earth in the 
Spherical harmonic mode Bar ates Spercleniis UNA eu s10 hl aeOn puma resc UreTl Cry, 


and the conductivity and is used to model observational 


data by comparison to the internal and external field ratios. 
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Lahiri and Price (1939) obtained an analytic 
solution for eq. 1.22 for a model whose conductivity 
varied as an arbitrary power of the radius. They fitted 
a five-parameter model to both Sq- and Dst-variations. Sq 
CONTAINS: INpontanmt Narnonics, at periods 6. os lic and 24 hy 
and Dst is a transient with useful amplitudes in the Fourier 
transform with periods up to several days. In the Lahiri- 
Price model the major features are the following: 

(1) The resistivity varies slowly from a value 
of 100 ohm.m near the earth's surface to about 5 ohm.m at 
a depth of about 400 km. 

(2) Between 400 and 800 km the conductivity 
increases rapidly with depth. 

Rikitake (1966) reporting work done around 1950, 
determined uniform core models which were compatible with 
available analyses of Sq, Dst, the geomagnetic bay (typical 
period 1 hour), solar flare effect (period of about 30 
min) and the 27-day period variation. He modelled the 
induction by means of a rapid fall of resistivity to about 
10 ohm.m at a depth of about 400 km. All the models 
obtained in the classical studies (Fig. 1.1a) show a steep 
rise in conductivity in the depth range 400 to 700 km, 
DUteitewasmnotepoOsSsiIDIe tO=Say COnnidentiy If, or where, 
the conductivity leveled off. Estimates of the deep 
conductivity have been made by studying the transmission 


through the mantle of the secular variation of the main 
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Figure 1.la Examples of conductivity models of the Earth 
obtained in classical studies (McDonald, 1957; 
Rikitake, 1960). The Price-McDonald curve 
published by McDonald (1957) combines his 
estimate of the lower mantle conductivity 
with one of the Lahiri and Price (1939) models 


Figure 1.1b The models of Banks (1969, 1972) and Parker 
(1970) for the radial conductivity of the 
Baer 
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geomagnetic field, which is generated in the outer core. 
Assuming that the cut-off frequency of the secular 


' as determined by Currie (1968)) is 


Variation (0725) yr. 
HoOtieansintrins ie property. but duerto. the fi tering 
properties of the mantle, a number of investigators (e.g. 
McDonald, 1957; Smylie, 1965) estimated the conductivity 
of the lower mantle to be 100 S.m7!, 

It was only about 1960 with the advent of high- 
speed computers capable of handling large data sets, that 
detailed analysis of the conductivity distribution in the 
earth became possible. Outstanding among the recent work 
is: Pthat of eBanks (1.969),. 19/72) land: Parker (19/08 1972). 
Banks examined the nature of the geomagnetic spectrum, 
the peaks as well as the continuum between them and found 
that the single spherical harmonic pe described the varia- 
tions quite well. An estimate of the response of the 
earth in the Br, mode over the frequency range 0.1 to .25 
day7! was obtained. A multi-layered model was fitted to 
these response data and produced a conductivity estimate 
aceeshownein eigen). bb aCBanks 4 61969) 8 Also shown im thts 
figure is an improved model obtained by taking into 
qccount caubetter fit ofuthe phase relations (Banks, 1972). 

Parker (1970) applied the Backus-Gilbert linear 
inverse technique (Backus and Gilbert, 1967, 1968), to 


the geomagnetic induction problem. He used as observables 
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(from the Banks (1969) data) the moduli of the geomagnetic 
response of the earth in the Pe mode at 35 frequencies 
between .01 and .2 taal The final model fitted to the 
data 1s) shown in Fig. tolb. “Parker s@model differs from 
Lidte Of banks sine tnet Ete gives a conductivity for the top 
400 km of the earth about one order of magnitude higher. 
The cause of this discrepancy could be the fact that 
Banks' data contain negative (nonphysical) phases at 
pernodss less than 5 dayse Of that Parker did not incor= 
porate the phase information in his inversion, or an 
inherent ambiguity of the modelling metnod Banks used, or 
more than one of these (Banks, 1972; Jady, 1974a,b). 

Bailey (1970) and Weidelt (1972) described 
exact inverse methods, which unfortunately require as 
Input precise data over a broad frequency range. Data of 
that quality are not yet available. 

invall studies discussed ane tits section, 
spherical symmetry has to be assumed to make the 
theoretical problem tractable. The short period limit 
for these studies is determined by induction effects 
caused by lateral variations in electrical conductivity 
ins the crust and Upper mantle. “In practise this means 
that only data with a period longer than about five days 


could be used in these studies. The outcome of this 


limitation at the high frequency end of the spectrum is 
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that the conductivity structure in the upper 400 km of the 
earth is poorly resolved (Rikitake, 1973). Studies of the 
top rew hundred Kviometres of the earth is) thus. best left 
to other methods like the direct current sounding methods, 
maqnetotelluric soundings and short-period geomagnetic 
deep sounding techniques not depending on the radial 
Symmetry requirement. 

Since the resistivity decreases to values less 
than 1 ohm.m in the depth range 400-800 km (Fig. 1.1b) 
even storm-time fields with periods of several days have 
highly reduced amplitudes at a depth of 1000 km. This 
limits geomagnetic induction studies by means of external 
fields effectively to depths less than 1000 km. The 
power of the geomagnetic variation spectrum indicates a 


3 = due to the 


steep rise at frequencies less than 10 “day 
geomagnetic secular variation (Rikitake, 1973; Banks, 
1969) which originates entirely within the earth. Even 

if very low-frequency external fields were available, the 
secular variation fields would mask them and would make it 
very difficult to use such very low-frequency external 


fields for study of the lower mantle. 


1.2.2 Direct-current sounding methods 
The two most widely used direct-current sounding 
methods as far as crustal-scale investigations are concerned 


are the Schlumberger and dipole techniques. 
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The Schlumberger electrode array consists of 
four colinear electrodes (Kunetz, 1966). To make a 
sounding the two outer electrodes, the contacts for 
current emission, are moved progressively away from the 
centre of the spread. The inner two electrodes, arranged 
Symmetrically with respect to the centre of the array 
measure the voltage drop due to the current flowing through 
the ground between the current electrodes. In the 
Schlumberger array the distance petween the current 
electrodes is always more than five times that between 
the potential electrodes and in deep crustal to upper 
mantle studies it can exceed 200 times the potential 
electrode distance. The ratio of the voltage drop to 
potential electrode separation thus approximates the 
electric field. 

In the dipole array the four electrodes are 
arranged in two pairs, one pair supplying current and 
the other measuring a potential difference, and the 
relative arrangement of the pairs determines the type of 
dipole array (Al'pin et al., 1966). The current electrode 
Haire is usually. fixed Tn) location, whilesia component) of 
the etectric field is mapbed as >a function of cistance 
from this current source with the potential-difference 
pair of electrodes. Normally the electric field is mapped 


away from the dipole source along one of the principal 
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directions. In an equatorial dipole sounding, the 
component of electric Tield parallel to the source axis 
is mapped outwards on the equatorial axis. For a 
horizontally layered earth - but only for this 
configuration - the resulting apparent resistivity curve 
as a function of separation is the same as one obtained 
from a Schlumberger sounding. When the electric field 
is measured at locations along the axis of the source the 
sounding is called a polar dipole sounding. 

In the direct-current sounding method the 


apparent resistivity Ps in onm.m is given by 


ae cky 
2 ee (2s 


With AY thes potential drop in Volts, 1 the current in 
Ampere and K a constant dependent on the electrode con- 


figuration and qiven by 


] 


where A and B denotes current electrode positions and 
Mand N potential electrode positions. The distances AM, 
AN, BM and BN are measured in metres. 

The classical method of interpreting dc soundings 
is by comparison of a field curve with a precalculated 


cataloque of reference curves (e.g. Orellana and Mooney, 
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1966). The curve-matching technique is, however, rapidly 
being replaced by computer-assisted interpretation tech- 
niques which permit more reliable results to be obtained 
with less effort (Kunetz and Rocroi, 1970; Inman, Ryu and 


Ward, 1973). 


1.2.3 The magnetotelluric method 


The maqnetotelluric sounding method makes use 
of the natural electromagnetic field of the earth as 
power source. The technique requires that simultaneous 
observations be made of orthogonal electric and magnetic 
variations in the electromagnetic field. Cagniard (1953) 
and Keller and Frischknecht (1966) provide excellent 
introductions to the theory of magnetotelluric fields over 
a plane-layered earth. 

With a coordinate system at the earth's surface 
having the axes aligned x north, y east and z down we have 


for nlane waves generated by sources in the ionosphere 


S 
Z = a (ie2an 


where Z is the characteristic impedance in ohms, EY the 
electracefield intensity (north) 1n)V/m and Hy is the 
magnetic field intensity (east) in A/m. In this section 
Capital letters Signi iy frequency-domain representations 
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In a homogeneous and isotropic earth the true 
resistivity is related to the characteristic impedance 


through the relation: 


(le2i5 3) 


where op is the resistivity in ohm.m, T the period in 
seconds and uw the magnetic permeability of the medium. 
In the case of a horizontally layered earth (one- 
dimensional earth) eq. 1.25 becomes that for the 


apparent resistivity 
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which is frequency dependent. 
In the presence of lateral changes in the 
electrical properties of the earth, however, each electric 
field component couples to both magnetic field components 


through a relationship of the form: 
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is called the impedance tensor (Word, Smith and Bostick, 


1970; Hermance, 1973). 


1.2.4 Geomagnetic deep sounding: linear array method and 
large two-dimensional arrays 


Classical geomagnetic deep sounding (G.D.S.) 
methods were developed mainly by Schmucker (1964, 1970a). 
In G.D.S. studies use is made of observations from a small 
number of three-component variometers, usually less than 
ten, recording simultaneously. Since the best arrangement 
of such a small number of instruments is often along a 
Straight line, G.D.S. methods are usually associated with 
small linear arrays. In regions removed from local 
anomalies the smoothness of the observed field indicates 
that the electrical conductivity varies approximately 
with depth alone, and is thus one-dimensional. In this 
Case it 1S the practise to compute the ratio Z/H of vertical 
to horizontal components for an incident field of given 
period T and spatial wave number k, at the surface of a 
layered conductive half-space (Schmucker, 1970a). The 
computed Z/H ratio is compared with the observed ratio for 
several periods at which data are available. These models 
are highly non-unique and the estimation of k is a major 
Gurrrculty,. [his problem Will De discussed 1) Chapter =o. 

In regions with lateral conductivity changes, 


the station-to-station changes in amplitude and phase of 
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the field components indicate the presence and configura- 
Eilon woTrsinterhalweonductive structure. “To omapy the conduc- 
Give structure’ through the use of different variation 
events, several events are recorded with the magnetometers 
in each of a number of positions successively. A measure 
of the response of the earth to the incident field is 
required, normalized to that field. Parkinson (1959, 1962) 
and Wiese (1962) took a first step with their arrow 
representations of response. Schmucker (1970a) introduced 
the use of a matrix of transfer functions, each of which 
gives the response of the earth in one component of the 
anomalous field to one component of the normal field. 

This method preserves the phase information. The esti- 
mation of transfer functions, however, requires computation 
of powers and cross-products between normal and anomalous 
field components. Unique separation of normal and anomalous 
fields is commonly impossible with a small array, and 
usually the field at one station is arbitrarily designated 
as representing the "normal" field. 

Large two-dimensional arrays of magnetometers 
came into use as a result of the development of an in- 
expensive portable three-component magnetometer by Gough 
and Reitzel (1967). -An advantage of an array is that it 
allows mapping of a single variation event over the whole 


array which may cover ice ie km?. iis peeves impor tant 
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when dealing with specially polarization-sensi tive 
anomalies. A Fourier transform refers to a precise 
period and sums the energy of the whole event at that 
bertod and for this reason the Fourier transform 
amplitudes and phases of variation components have 
proven good parameters to map (Reitzel et al., 1970). 
The data analyses involved in array studies will be 
discussed in Chapter 2. 

Magnetic fields observed over part of the 
earth may be separable into external and internal parts 
by use of surface integrals. Porath, Oldenburg and Gough 
(1970) used surface integrals to separate fields obtained 
with a two-dimensional magnetometer array. The separation 
can be done for instantaneous fields in the time domain 
or for a given frequency in the frequency domain. 

Separation OT Variation "fields over a limited 
portion of the earth, however, has some short-comings. 
Fields whose dimensions are larger than the dimensions of 
the array are not separable by the method. In the middle 
latitudes the source fields have dimensions that are larger 
than the array dimensions and cannot be separated. These 
limitations led Gough (1973b) to state: “The separation 
exercise described by Porath et al. (1970) was valuable, 
but it may be doubted whether it should be repeated." 

Quantitative interpretation of conductivity 
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i “7 ines 
f iy } 
; 2 
' 
* 
4 
+ Sttrit. 
Ld in ( 
if 
! ! i i 
' i ft 
: 1988 “7en, OARS sry. tm 
A 
Ll it : ; 
iA ote 
a +i dha Ley eit ee) 


as en: aa rrial? # 
eat ys vevens oA RETR 
daa oth wednsy ‘hill ich oo i 
tor IP pte, patty bah 2 va | Pome 


| ied id -shenethe boas ave 
" we _ 


ia Th? tae ys 
ane 
t ” 
’ i 
ris Fi 
Cte © 1b brik ae 
i ou afi 
ae | 
Poh 7 
4 1m] < | f { 
reg ‘u 1 eR ane 


BIDS ENE rt) 


AGhaniead Saslaiworee uate ce a4 


232 Ar wo i 


Kings nits +8 | pha 7 


“hea 45°F) Vo) 2iyata bs "a Fou 


roman 
te ih bl Fey ae ov gp sit: 
| ce ons. abe Pe. 


4< we 7G4 as sah ab 32 


a 


ca 


to the anomalous fields, normalized with respect to the 
normal (inducing) fields. Analytic methods of modelling 
deal only with conductors of a few simple geometrical 
Shapes. d'Erceville and Kunetz (1962) obtained exact 
solutions for a semi-infinite conductor with a vertical 
plane of discontinuity. Rankin (1962) solved the problem 
of a dyke model exactly. Weaver (1963) and Weaver and 
Thomson (1972) produced analytic approximate solutions 
for a non-uniform earth with an overhead line current. 

Analogue scaled models have been used to study 
conductivity structure (Dosso, 1966; Dosso and Jacobs, 
1968; Hermance, 1968). 

The most successful approach to the problem was 
through numerical methods. In two-dimensional cases it is 
useful to consider anomalous fields which are E-polarized 
or H-polarized, In the E-polarized field, the induced 
currents flow parallel to the infinitely long dimension 
of the structure. The induced field is also E-polarized. 
The induced currents of the H-polarized field flow at 
rignt angles to the long dimension of the structure and 
their magnetic fields are not observable at the earth's 
surface. 

Of the numerical methods used for modelling of 
arbitrary shapes three have been widely used, viz. the 


transmission line analogy method (Madden and Swift, 1969; 
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Wright, 1969), the finite difference method (Jones and 
Priee, 1970;— Jones and Pascoe, 1971: Pascoe and Jones, 
1972) and the finite element method (Coggon, 1971; Reddy 
and Ranking 19/3). 

Geomagnetic deep sounding techniques proved to 
be especially well suited to detect lateral changes in 
conductivity in the crust and uppermost mantle (Porath, 


TIITEM GOUGH 919 74)2 


lec sue CONCUCtTIVAItYRdiStrabutioneain: tie «crust and Unpermos & 
Many attempts have been made in recent years to 
determine quantitative electrical resistivity profiles for 
the crust and unpermost mantle by means of surface measure- 
ments (e.g. Keller, 1971). The vast majority of these 
investigations were conducted by studying variations in 
the natural electromagnetic field of the earth, the 
magnetotelluric method being the most widely used technique. 
It is widely recognised that the stable continen- 
tal crust can be divided into three broad units as far as 
Htsselectrical conductivity ws concerned (e.g. Keller et ale, 
1966; Dvorak, 1975). The uppermost of these units is the 
conductive surface layer which consists of sedimentary or 
decomposed igneous and metamorphic rocks sometimes containing 
water. The underlying second unit is highly resistant (of 
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the order of 10° ohm.m in stable continental regions (Keller 
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et al., 1966) and is made up of fresh igneous and meta- 
morphic rocks which although buried, are still at low 
temperatures but extend to depths where pressures are 

high enough to close cracks and exclude water. As the 
temperature increases with depth, the resistivity decreases 
and the third unit, which is distinctly more conductive 
than the second, is reached. This unit has a resistivity 
of about? ,000 ohm.m for the stable continental model 
(Dvorak, 1975). Another interesting feature in the 
conductivity profile is that there is no decrease in 
resistivity associated with the Mohorovici€ discontinuity 
(Dvorak, LOTS ie tease Worthy. OT note that wes s CIV 1Cy 
ranges through at least 6 orders of magnitude in the earth, 
whereas seismic wave velocities, for instance, do not cover 
aaractor of 10" 

In a conductor-resistor-conductor sequence of 
units as in an idealized crust, the magnetotelluric method 
is sensitive mainly to the thickness of the resistant zone 
and gives a poor indication of its resistivity (Madden, 
1971). Direct current sounding methods which use a 
controlled artificial source of electrical energy, however, 
are sensitive to the transverse resistance Ty (thickness- 
resistivity product) of the resistant zone, especially if 
the surface layer has a large longitudinal conductance S, 


(thickness-resistivity ratio). If, however, areas can be 
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found where the conductive first zone is absent or where Sy 
is small it should be possible to determine the true 
resistivity and thus the thickness of the resistant zone 

or at least to approach the true value. The direct current- 
resistivity methods are thus from a theoretical point of 
view the most suitable methods to study the resistivity 
distribution in the resistant portion of the crust (Madden, 
opis) 

The dc and MT methods are probably comparable as 
regards the investigation of the third unit because the 
response of both methods depends on the longitudinal 
conductance S 3. 

Both the dc and MT methods should ideally be 
used in areas where the conductivity distribution 
approximates a horizontally stratified case. Although 
local geomagnetic deep sounding studies can be, and have 
been, used for studying horizontally stratified situations 
(Schmucker, 1970a> Poratn et al... 19/713 Kuckes, 19734, 
Camfield and Gough, 1975) it is more suitable for mapping 


lateral conductivity changes. 


1.2.6 Results of some magnetometer array studies 


Linear magnetometer array work discovered 
several important conductivity anomalies. Examples are: 


the Rio Grande anomaly near the Mexican border in the 
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Midwestern United States (Schmucker, 1964, 1970a), the 
extensively studied Japanese anomaly (Rikitake, 1966; 
Rikitake and Honkura, 1973) and the conductivity anomaly 
under the Peruvian Andes (Schmucker et al., 1967). 

The first two-dimensional array studies were 
carried out in western North America by the University of 
Texas at Dallas and the University of Alberta (Reitzel, 
Gough, Porath and Anderson, 1970; Porath, Oldenburg and 
Gough, 1970; Porath and Gough, 1971; Camfield, Gough and 
Porat, 9/1 =Porath, Gougn and Camfield. 197)3 sGamfield 
and Gough, 1975). These initial studies were followed by 
studies with arrays from either of the above institutions 
(e.g. Porath and Dziewonski, 1971la; Alabi, 1974; Alabi, 
Camfield and Gough, 1975). 

Large arrays were also operated in Australia 
(Gough, McElhinny and Lilley, 1974; Bennett and Lilley, 
T9737 eeScotlanducsrk eteale. 1974). Ethiopia Berk told 
et al., 1974), Kenya (Banks and Beamish, 1974) as well as 
South Africa and South-West Africa, Botswana and Rhodesia 
(the present study, Gough, De Beer and Van Zijl, 1973; 

De Beer et al., in press). 

Comprehensive reviews of magnetometer array 
studies have been published by Porath and pziewonski (1971b), 
Gough (1973a,b), Frazer (1974) and Lilley (1975). In this 


section I will therefore only summarize some of the more 
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important results from the above two-dimensional array 
Studies and indicate the correlations between the electrical 
conductivity anomalies and other geophysical anomalies. 
The complicated conductivity structure of 
western North America has been studied in great detail by 
the arrays that operated in that area. Gough (1974) gave 
a comprehensive review of the results. It was found that 
highly conductive mantle material which has to be at least 
IGOR Meth tees Gnea resistivity. of soronm. mm cunderiieswthe 
Basin and Range Province and ridges of still lower resis- 
tivity (or greater thickness for a given resistivity) 
underlie the Wasatch Fault Belt and the southern Rockies. 
The upper mantle under the Great Plains proved to be more 
resistive and may approximate a sub-shield type. The 
conductivity under the Colorado Plateau is of the Great 
Plains type or intermediate between this and that beneath 
the Basin and Range. Gough emphasizes in his review the 
ambiguity of interpretation of depths of structures 
associated with the models for the section through the 
two belts of high conductivity on either side of the 
Colorado Plateau. The array results could be satisfied 
equally well by undulations in the depth to the highly- 
conducting mantle (giving depths of between 100 and 400 km) 
OM Dyevaritations in thickness of “a sconducting layer inj the 


upper mantle separated from the deep mantle. The second 
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interoretation 1s more attractive because it allows a 
comparison with the low-velocity zone. The seismic velocity 
models for western and eastern North America show a 
striking difference in thickness and intensity between the 
two portions of the continent. In the west the depth 
extent of the zone corresponds well with the required 
thickness of conducting material, while the much 

smaller low-velocity zone in the east is compatible 
with the fact that conductivity anomalies in eastern North 
America appear to be controlled by crustal rather than 
mantle conditions (Garland, 1975). High heat flow in the 
western United States also shows a good correlation with 
tive tI Strroucion of structures swith high conductivity. 
Gough (1974) points out that there can be no serious 

doubt that the highly-conductive regions of western North 
America are conductive because they are hot. 

The North American Central Plains magnetic 
variation anomaly is one of the most striking elongated 
anomalies (Gough and Camfield, 1972; Alabi, 1974). The 
conductive body underlies an area where the basement rocks 
are poorly exposed and its presence was indicated for the 
first time by magnetometer array measurements. The above 
methods have shown that there is supporting evidence for a 
metamorphic belt (graphite schist belt) in the crust and 


faulting in the basement under the Palaeozoic sediments. 
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This anomaly, first detected near the Black Hills, is also 
of the current concentration type which means the anomalous 
fields are those of currents concentrated in a body of 
high conductivity which joins large undefined regions of 
the Earth ini wnich the induction occurs... This effect 
precludes two-dimensional modelling of the anomaly. 
The first array study in Africa was over the 
Afar Depression and Main Ethiopian Rift in Ethiopia (Berktold 
et al., 1974). Analysis of the spatial inhomogeneity: of the 
magnetic variations near the magnetic dip equator which 
crosses the southern part of the Afar Depression indicates 
a resistivity of less than 10 ohm.m at a denth of 100-300 km 
in this area. Higher electrical conductivities are also 
indicated under the Afar Depression than beneath the adjacent 
plateaus. South of 12°N, however, Parkinson arrows suggest 
higher conductivity beneath the southern part of the Western 
Plateau than beneath the Afar Depression. The above authors 
suggest that the low resistivities are temperature related. 
Banks and Ottey (1974) found a similar situation 
from a linear array study across the eastern branch of the 
East African Rift Valley near the Equator. They concluded 
that the type of model that fits the experimental data 
involves a strip of high conductivity material at a depth 
of no more than 20 km beneath the floor of the rift with, 
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moderate conductivity at a depth of 50 km beneath the 
eastern flank of the rift. The two regions of high 
conductivity are interpreted as zones of partial melting in 
the upper mantle. An associated low velocity, low density 
zone is indicated by the gravity and seismic experiments 
which support the interpretation of partial melting in 

the upper mantle. Additional strong evidence in favour of 
such an interpretation is the recent volcanic activity in 
and around the rift valley and the associated geothermal 
activity. A two-dimensional array was operated in the 
Same region during 1972 (Banks and Beamish, 1974). Data 
from the east-west line of instruments roughly along the 
equator appear to confirm the linear array results and 
indicate a two-dimensional conductive structure beneath 
and to the east of the Rift Valley. The results from this 
array also indicate a three-dimensional structure under- 


neath the area to the south-east of the 1971 line. 


1.3 Geological background to the investigation 


The tectonic pattern within the Precambrian 
continental basement in southern Africa is characterized 
by a network of intersecting mobile belts (Fig. 1.2). 
These mobile belts surround and separate from one another 
cratonic areas which have remained undeformed for at least 


the last 2,500 Myr (Anhaeusser et al., 1969). 
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Eburnian-Haubian, 


Orogeny 


Sketch man of shields older than 2,500 Myr (shaded 
areas) and mobile belts in southern Africa. The 
numbers from 1 to 8 indicate the centres of deep 
Schlumberger soundings listed in Table 1.1 and the 
heavy lines indicate telepnvhone and power lines 
used as current emission lines in these soundings. 
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Clifford (1972) argued on the basis of available 
geologic and geochronologic data that the major part of 
the African continent is a segment of primeval crust 
more than 3,000 Myr old, which has been affected by 
several major orogenies. At least two of these orogenic 
events are older than 2,500 Myr and the best evidence 
for their occurrence is preserved in the cratonic nuclei. 
The major younger orogenic events in southern Africa 
occurred as a sequence of polyepisodic orogenies 
1,850 + 250 Myr ago (Eburnian and Huabian orogenic 
episodes), 1,100 + 200 Myr ago (Kibaran orogeny), 

600 + 100 Myr ago (Damaran-Katangan or Pan-African orogeny) 
and during middle Palaeozoic - early Mesozoic time when 
the Cape Fold belt was formed (Clifford, 1974). 

The recognition, delineation and structural 
analysis of these tectonic units have been the topics of 
many studies. Outstanding in the early work in southern 
Africa is the interpretation of Holmes (1951) who combined 
tectonic studies and the then new study of geochronology to 
discern internal features within the Precambrian terrains. 
Similar studies by Cahen (1961), Nicolaysen (1962), Cahen 
and Snelling (1966) and Clifford (1968) followed. 

Detailed analyses on limited areas have been made by, among 
others, Nicolaysen and Burger (1965), Vail (1965), Clifford 
(1967) eiVotnee teal 1908.) . 041 isoppee tal 01969) Vain and 
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Snelling (1971), Van Breemen and Dodson (1972), Wilson 
andenarvisoMm Cloyd) and Clifford etal, (1975). 

These studies established the present knowledge 
about the distribution and interrelationship of the differ- 
ent SLUUG UU Laut Nts: dS  Cepicted: 1nmimiGnlsc. 

A =brvet description cof ithe structural units 
qs snown in rig. Isc will waciilitate tater discussion. 

The Kaapvaal and Rhodesian cratons consist largely of a 
granite and gneiss terrain incorporating early Precambrian 
greenstone belts. The whole suite is older than about 2,600 
Myr, the age of the last major granitic and pegmatitic 
intrusion (Nicolaysen and Burger, 1965; Martin, 1969). 

The ancient nuclei within these 2,600 Myr old terrains 
GONMSUS tsa. On CIM crus Ver Ss UMteS: Of GYanibtdc rocks: and Scns t 
belts of metasedimentary and volcanic rocks with ages up 

to 3,400 Myr (Van Niekerk and Burger, 1969; Allsopp et al., 
1969; Vail and Dodson, 1969). 

The Limpopo Belt consists of reworked Archean 
granite-greenstone terrain with an early Proterozoic cover 
sequence infolded and metamorpnosed with the basement (Mason, 
1973). Van Breemen and Dodson (1972) found that although 
the main folding and metamorphism dates back to 2,/00 Myr, 
a strong thermal event affected the Limpopo Belt about 
2,000 Myr ago. 

The Namaqualand-Natal Belt is constituted of 


granite gneiss, together with remnants of the metavolcanic 
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and metasedimentary rocks of the Kheis System that are 
older than 2,600 Myr. The near-equality of ages of 900 

to 1100 Myr yielded by a wide range of minerals is believed 
to reflect a period of intense reconstruction of older 
basement (Nicolaysen and Burger, 1965). 

The Zambezi Belt comprises a zone of high-grade 
metamorphism around the north-western and northern margins 
of Rhodesia. Down-faulted Karroo (Permo-Triassic) rocks 
in the Zambezi Rift Valley cover the metamorphic rocks 
of the belt for tne largest part. Geochronological results 
indicate an age of about 500 Myr for the latest thermal 
event in the belt (Vail and Snelling, 1971). In the west 
the Zambezi Belt disappears under the Cainozoic Kalahari 
cover and in the east it joins the north-south trending 
Mocambique Belt. The Mocgambique Belt consists of high-grade 
metamorphic rocks with associated granites and peqmatites. 
This belt is also characterized by ages between 400 and 
700 Myr (Clifford, 1967). The Damara geosynclinal belt 
consists of a miogeosynclinal facies in which dolomitic 
limestone predominates and an eugeosynclinal facies con- 
Paina mMatie Volcan cs. and: consisting Of predominantly 
pelitic sediments. The most recent orogenic event in the 
eugeosyncline is indicated by consistent radiometric ages 
in the general range 450-550 Myr (Clifford, 1967). The 


similarity in the age patterns between the Damara geosyn- 
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cline, Zambezi Belt and Mocambique Belt suggests a tectono- 
thermal kKinshipeot the? zones. The exact structural relation- 
ship between the Damara geosyncline and the Zambezi Belt 

1s, however, obscured by the Kalahari sediments in Botswana, 
southern Anggla and western Zambia. 

The youngest of the major tectonic events 
observable in South Africa is represented by the Cape Fold 
Belt. This fold belt consists of the less severe north- 
south trending Cedarberg folds in the west and the main 
east-west trending fold system in the south. The folded 
sedimentary strata range in age from Late Ordovician 
(Cramer et al., 1974) to Late Permian (Haughton, 1969), but 
theamainypnases of Tolding occurred in iyiassic times. 
Several authors suggested a continental collision origin 
rotathestonmatwoneot= this fold=belt (De Swardteet alk, 1974; 
DemBeersetual..9 1974 Martinis 19745 Rhodes.) 1974)% 

Clifford (1972) concluded that these mobile zones 
are floored by older crystalline basement and that no geo- 
‘Synclinal sequences were deposited on extensive areas of 
oceanic crust. In large segments of almost all these 
orogenic zones, the orogenic event is entirely recorded as 
rejuvenated floor rocks. Shackleton (1973) makes the 
iMDO MEAN LPO te tilatmthereraresno large offsets of earlier 
structures across these belts. This implies that where two 


Momllewpeltsecrosss there Nias been Jittle, if anya crustal 
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shortening in the younger orogenic belt and thus little, 
if any, relative motions of the stable cratons on either 
Side of this’ belt, 

Despite all the geochronologic-tectonic studies, 
the information regarding both extent and age of the units 
oes il ievambiquous,, (ihe difficulty of Interpreting isotopic 
age measurements as dating either tectono-thermal events, 
or merely regional cooling through uplift without orogenesis, 
inhibits the correlation and paeaone ean of many of the 
active pelts (Vail a 1900).8 The overprinting of isotopic 
events on stable cratonic blocks or earlier orogenic 
belts without the apparent accompanying metamorphism may 
also mask the aerial extent of the units. Undisturbed 
younger sedimentary cover further obscures vast areas of 
the Precambrian—of Africa, especially since the cratons 
have provided basins for considerable sedimentation during 
Upper Palaeozoic and Mesozoic times (Kennedy, 1965; 
Haughton, 1969). The Palaeozoic-Mesozoic Karroo strata 
obscure in this way the transition from the Kaapvaal 
shield to the Namaqualand-Natal Belt. The Cainozoic 
Kalahari beds in central soutnern Africa similarly impede 
a definitive structural correlation between the Damara 


geosynclinal belt and the Zambezi Belt. 
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1.4 Results of ultra deep Schlumberger soundings in 


southern Africa 


During 1967 the National Physical Research 
Paboral omy (Nera tetiOh the eS leRee Sout HAY Ca 
embarked on a programme of ultra deep Schlumberger 
Soundings to investigate the electrical conductivity 
Structure of the crust and uppermost mantle in southern 
RErical Van. 21 jl ewe Goose Valet} leet ole. 1O70Re Van 215 1 
and Joubert, 1975). The positions of the eight soundings 
carried out since the start of the programme are indicated 
Vie ae lace aD ve” ie ledi ves. thesstectonic, oroav ince on 
which the centre of each of the soundings was situated, 
the maximum current electrode spacing and the transverse 
resistance of the resistive part of the crust. In these 
Soundings telephone lines and power lines were used as 
current emission lines. 

The programme proved to be very successful in 
establishing a resistivity profile for the major 
SUructurale UnIts in the country. lable 1 shows: that 
the transverse resistance of the resistant part of the 
crust in the Kaapvaal and Rhodesian cratons is from a 
half to two-thirds that in the Namaqualand-Natal Belt. 
This is most nrobably due to a more complex metamorphic 


histowvernethe teattem reqion (Van Zigl et al, 1970). It 
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Tectonic provinces on which centres of Schlumberger soundings 
were situated, maximum current electrode spacing and trans- 
verse resistance of second unit 


Maximum 
current Transverse 
electrode Resistance of 
Sounding No. Tectonic spacing second unit, Te 
(elpdesle cs) Province (km) (ohm. m? ) 
] Namaqualand- (AL) ial ees lee 
Natal Belt 
ra Namaqualand- 330 Sie ox 10° 
Natal Belt 
3 Namaqualand- 139 not determined 
Natal Belt 
: 8 
4 Kaapvaal 600 5 x 10 
craton 
8 
5 Damara 300 1x 310 
Eugeosyncline 
; 8 
6 Rhodesian 450 bo <0 
Cia con 
8 
7 Kaapvaa|l 450 Oo en 0 
eraton 
8 Mocambique hoo not determined 
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also shows that the transverse resistance of the resistant 
unit in the Damara eugeosyncline is in turn about five 
timessiess than what of the cratons., —Fidure 1.3 depicts 
three of the latest sounding curves on the cratons and 

the Namaqualand-Natal Belt together with the best fitting 
models. The conductor-resistor-conductor structure for 
the stable continental crust is very clear. The models 
Show that the crust has a maximum resistivity not 


2 ohm.m at a depth of less than 19 km. 


exceeding 10 
Below this highly resistive layer the resistivity then 
decreases progressively to reach a value of 700-1309 ohm.m 
at a depth of 16-20 km, after which the resistivity does 

not decrease appreciably until a depth of at least 200 km 

is reached. No obvious change in resistivity thus exists 

at depths comparable to the Mohorovitié discontinuity 

which is at a depth of 35 to 37 km in the Kaapvaal shield 
(Green and Bloch, 1969). 

Despite the good results obtained, the deep 
Schlumberger sounding technique has the big disadvantage 
that the availability of long emission lines dictates the 
study area. JThe investigator thus has very little control 
over the siting of a sounding with respect to the geological 


province, and can only make studies where the necessary long 


lines can be obtained. 
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Experimental Schlumberger sounding data (discret 
symbols) and curves for best fitting models for 
soundings on the Namaqualand-Natal Belt (curve 1 
Rhodesian craton (curve 2) and Kaapvaal craton 
(curve 3) are shown in the upper diagram. These 
soundinas correspond respectively to soundings 
numbered 2, 6 and 7 in Table 1.1. The lower 
figure depicts the best fitting resistivity 
models for the three cases. The thin weathered 
surface layer is omitted. 
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In order to expand the information obtained by 
the Schlumberger soundings the N.P.R.L. thus introduced 
a programme of magnetometer array studies to investigate 
lateral changes in conductivity in the crust and upper 
mantle in southern Africa. The analysis of the magnetometer 
array data and the interpretation of the results obtained 
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CHAPTER. t 
DATA ANALYSIS AND INTERPRETATION TECHNIQUES 


201" lhewid ata 


The analysis of magnetometer array data has 
been described in great detail by Anderson (1970) and 
Camfield (1973) and has been discussed by Gough (1973b). 

The Gough-Reitzel magnetometers record the H 
(magnetic north), D (magnetic east) and Z (vertical down) 
Components, inthe period range 2 min < 7 =< 1-day, -on 
35 mm film. These films are examined and sections of 
record showing clear magnetic disturbance events are 
chosen for analysis. In the present study the sections of 
35 mm film recordings of the selected events were printed 
electrostatically at ten times magnification, edited and 
all deflections with a period of less than two minutes 
smoothed by hand. The smoothing was Carried out to prevent 
alvasing during the commercial digitizing at one minute 
intervals from the prints. Full-scale plots of the 
diaqital values were superposed on the prints and carefully 
checked to ensure that the digital values were within one 
tracewidth from the recorded curve. Values with larger 
deviations were corrected and in some cases the entire 
record was rediaqitized. [he digitized trace ordinates 


were converted to nanotesla (gammas) using the scale factors 
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derived from the field calibration of the magnetometers. 

A first order correction for a small interaction between 
the H and Z magnets in the magnetometer was also made at 
this stage (Camfield, 1973). The H and D (local magnetic) 
co-ordinates were also rotated to X and Y (geographic) 
co-ordinates. 

With the data in this form stacks of magnetograms 
were prepared and the first qualitative interpretation of 
anomalous variation fields carried out. This type of 
interpretation will be discussed in Chapters 3 and 4 for 


actual array studies. 


PwcCammOpeCLdileaial VSise and maps aOt srOUnietes DEC tral 
components 


Quantitative analysis of magnetic variation data 
is best performed in the frequency domain, because the 
depth of penetration of the magnetic field over a given 
structure is controlled by its frequency and spatial 
wavenumber. In this work simple magnetic disturbance 
events such as geomagnetic bays are considered as complete 
transient events with definite onset and finish times. 
Phethat 1s the case Fourier transformation otf the transient 
gives the complete and exact representation of the set of 
infinite waves of which the event is composed. The spectral 
density function at any frequency can then be considered to 
reveal the response of the system at that frequency. This 
allows examination of simultaneous frequency responses over 


a whole array area. 
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The definition of the Fourier transform used is 


eae: 
F(w) = fit yee Ode 
-T/2 
ee: 
= j Pte Onde 
-T/2 
wi th ie te Sol nc (See e.g. 


Tomcmen eee 


In the digital case we have 
Xa Orie Crmcemet se Nile Wet Meets =) pte ty 
Therefore 
Pen IKAG 
2N-1 ~ 2NAt 
ay = At 2 X. e 
je0) 
2N-1 ' 
_ : pee Tew 
= At a X{cos NAT 


The actual Fourier analyses 


performed with the aid of the IBM SSP 
This subroutine calculates cosine and 


and By for an input Xj of data points 
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where AMS UG lieeie eg ONAL 
2N-1 
ad 1 jk 
AL ait nA Se cos 1 
J (gran) 
2N-1 
and By = v r X. sin THK 
eth 


(See e.g. Jenkins and Watts, 1968, p. 19). 


If the real part of FY inv eq..-2.2- 1:5 denoted 


as Ey 1% is clear that 


ates AL > 2NAt 
k Z 
(255) 
nea 
2 
and similarly for the imaginary part 
F Bat. 
je ee gtk 
By = (20561) 


To relate the coefficients of RHARM to the Fourier 
trans torm as derined by eqs = 2.1 and 2.2 tne Ay and By of 
@qso2.4 are thus multiplied by t/2 (Bennett, 19/72). 

The Cooley-Tukey (1965) fast Fourier transform 
algorithm on which the RHARM subroutine 1S based requires a 
work vector of 2” points (M integer) for maximum computing 
efficiency (Gentleman and Sande, 1966) and for this study 
he was taken as 2048. The data lengths analysed were much 
shorter than this, but the remaining points of the work 


vector were filled with zeros. Oldenburg (1969) showed 
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that the addition of zeros in the work vector has no effect 
on the frequency spectra. A base line and a linear trend 
were removed from every event analysed and a sei taper 
varying from.zero to one was applied to the first and last 
ten minutes of data to remove fictitious high frequency 
components due to abrupt truncations at the beginning and 
end of the event. In the cases where the entire event was 
Fourier analysed, this taper was probably an unnecessary 
precaution, because every event was selected so as to include 
a smooth interval before and after it. 

The Fourier amplitudes as used for contouring 
are calculated as the square root of the sum of the real 
(cosine) term squared and imaginary (sine) term squared. 
These amplitudes are quoted in nT/c/ min since the digitizing 
interval At was in all cases one minute. The phase is the 
arctangent of the sine term divided by the cosine term. 
As pointed out by Bennett and Lilley (1972) this gives phase 


lag since if the input is a frequency component cos (wt-$¢) 


then 
COS toon = COS tECOS OmeeS Inputs (ed 
ve RO UNMetee stale Ce lel e SHEN ON a 2.7 
givinguarctaw Aronmer cosinenterm °° °@" \cosng! - 9 (2-7) 
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A representation of a variation event at a 
selected neriod requires six maps, which may show 
amplitude and phase for the three orthogonal components 
Keel 20% COSTING -andasine Coetti cients: tor tie tnree 
components. These maps of Fourier transform amplitudes 
and phases were introduced in the first array study by 
the Alberta and Dallas groups (Reitzel et al., 1970; 
Gough, 1973b) and serve two main purposes. Firstly they 
give a great deal of qualitative and some first order 
quantitative information and secondly form the basis of 
any further quantitative interpretation. 

A very important point is that contour maps of 
Fourier spectral components show both source and response 
functions and in this way differ drastically from more 
familiar maps such as aravity anomaly maps. If local 
anomalies apnear consistently with approximately unchanged 
positions for different variation events, one can with 
some confidence assign them to internal conductive structure, 
unless the array is under either the equatorial or auroral 
electrojet. 

The use and value of these contour maps will 


be demonstrated in Chapters 3 and 4. 
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The sources of error in the spectral components 


which are used in constructing the contour maps and which 
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alsostorm the inouteto all other interpretation in the 


frequency domain are the following. 


ZasubseGalibration errors 


At least two and sometimes three estimates of 
the calibration factors were obtained for each magnetic 
component from the calibration marks at the start of 
every film. Two estimates were used in the cases where 
the first two factors obtained differed by tess than 17: 
Porrthe: Tinst array 432 "and for tielisecond /6270f the 
CaltDrations fell “into this croup. Lf the tfirs i -two dif fered 
by more than 1% a third estimate was obtained. The 
differences were, however, always less than 3% which can 
be taken as a 95% confidence limit or two standard 
deviations. An estimate of the standard deviation of tne 
Calibration tactors vs thus 1.52 with corresponding errors 


in the amplitudes of spectral values. 


eee al Mim ernn6 3S 


With the aid of the time calibration marks at 
the beginning and end of every film the drift rates of 
the Accutron timers were determined. This gave relative 
times between instruments to within a minute. An additional 
CHeckeonm the timing. and linearity of tne drift rate was 
made within every selected event by examining Pi2 micro- 
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minute) over the whole array (Rostoker, 1972). The 
relative starting times of the magnetic disturbance 
events are thus known to within a minute and in most 
cases within 30 seconds. 

The errors in spectral phases due to timing are 


given by 
Ad = fAt + 360 degrees (2078) 


where At is timing error (taken as timing interval, one 
minute) and f = frequency (cycles/min). With the timing 
error one minute the corresponding phase errors are thus 
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2N-1 _injk 
es N 
E{d, | = | » d, e 
j=0 
ON-1 sh 
r ELd;] e 
j=0 
= G) G29) 
é 2 * 
and Variance id) J = E{d dy] 
= of (2.10) 


where dy is the complex conjugate of dy (Otnes and 

EMOCH GON malG 7 com Dee eo) ea ii Suwa bel ellSomgive «the vdaiudital 
errors in individual sine and cosine coefficients. The 
corresponding errors in amplitude and phase at a specific 
frequency can be calculated from these errors by normal 
propagation of errors formulae (e.g. Bevington, 1969). 

An estimate for a, can be obtained in the 
following way. As discussed before, the digital values 
were determined to within one trace width (about 2 mm) 
on the xerox enlargements. This means that we know 
with 95% confidence that |d | <2 mn or more specifically 
that om corresponds to 1 inn On the xerox prints. nis 
implies for average sensitivities of 1.001 nT/mm_ for X, 
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present array studies, that the standard deviation a 1s 
PCO en whore we Cole nil tomy and 2/19omth fom 2.  erom 
eq. ec. 101 tollows that for a four hour long event the 
Standard deviations of the cosine and sine coefficients 
eneCan oOn Mun Cy lt et Ole Age lO On ta ce MithietO leased. lilac 
ee Cyan 0 YoeeLec 
Since the longer periods usually contain more 

DOWeY. Chis constant standamd deviation for all meriods 
Will Mave a larder effect on the accuracy of the short- 
peGiod  FOUrTer spectral valves. Slt 1s Inrormative to 
consider examples from the next chapter. Figure 3.8 
Shows Fourier transform amplitudes and phases at period 
CAMMAMU Ce Sec Olomarcs UDSHLOMMasciats OCCUrre ds lon 30) = 22.30 “Galle: Ie 
on Auqust 31, 1971. The amplitudes are given in 
WS aig os 1! and the pnases in minutes. The average 
Standard deviations as percentage of the X, Y, and Z 
amonlitudes are respectively 12%, 24% and 25% and 3%, 4% 
and 6 fOr the pases. | Frqure 3.10 Shows Fourier trans torn 
amplitudes .tnen t/-c/min x 192 ) and pnases (in minutes) at 
aepeGlodrof 20 Minx <againetOr Une event Of AugUS Tag IE 
In this case the average standard deviations in X, Y and 
Lamplitudes due to digitizing errors are respectively 142, 
ma tGune cama (Gaal men ycnraNidiec 26 O Gites Vince olas ese 

ise clearetnat ot Shomt periodss ties digi ti zimg 
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and phase, while at lonager periods the magnitude of the 
dvaTenzind errors 1S comparable to that of the errors. in 


calibration and timing. 


Pea POlavizavioneore tne. norizontal, flea 


ies pelariezacionvof thertioruzontall afield. 1s 
important in geomagnetic studies, especially when induction 
in three-dimensional bodies is involved. With the aid of 
monochromatic wave theory the polarization of the horizontal 
field can be determined for any frequency component from 
the frequency spectrum of a transient. The horizontal 
field polarization is highly relevant to the interpretation 
of Fourier component contour maps for regions of 
complicated conductive structure (Gough et al., 1974). 

The polarization parameters for the horizontal 
variation fields were calculated following the monochromatic 
wave theory discussed by Born and Wolf (1959, p. 34). Ihe 
adantation of the formulae for the geomagnetic case (the 
X-Y axes of geomagnetism are reversed from the conventional 
2-D x-y axes) is described by Lilley and Bennett (1972). 


If a disturbance has horizontal components 
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then the point (X,Y) traces out an ellipse with the 


passaae of time. 


The anqle a, (0 < a < ©) is defined such that 
ia 
oD 
tan a = — (272) 
= 


The principal semi-axes a and b of the ellipse 
and the angle wv (0 < wy < wt), the angle measured clockwise 
between the major axis of the ellipse and the geographic 


east are specified by the formulae 


- aE pe = a‘ 1 eG 
tan 2us=) (tan 2a) cos 6 (Zeal) 
Sineeye= Stn co} sim 6 

where OR Un iCmenl fh BCOS (Gr 20, 


ic eon mn @C OST oma), 
a= Oba feecos; 6 = 0% cand a, > aos 
Ve= nse. th ecos 6 = .0) “and a, < a5, 


and y, (-m/4 < y « 0/4) 1s an auxiliary angle which 


determines the shape and orientation of the ellipse, 


tan xy = + 2 CMD) 
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Ltismine =< 05 the ellipse 1s ‘described in a 
clockwise sense and if sin 6 > 0, the ellipse is described 
anticlockwise. With sin o-= 0. or 4 (e7 ther a, or ap eS 
zero, the ellipse degenerates to a straight line. The 
values ays ay and 6 as functions of the angular frequency 
w come from the Fourier transforms of tne selected events. 

IMetnis Study Tt was found thar the horizontal 
field polarization varied somewhat across an array area, 
more at the sharter periods tnan at the longer periods. 
This is most probably caused by changes in the amplitudes 
and phases of the horizontal fields due to different 
NaUCTION eGftects across. the array. Another contributing 
factor may be the errors in tne Fourier spectral 
components. The horizontal polarization ellipses shown 
on the contour diagrams in Chapters 3 and 4 are ellipses 
from stations removed from the main anomalous structures. 
These ellioses wil? thus show the polarization in the 


more “normal” regions of the array. 


2. oo lrans fer functions  andocinauction Vectors 


Transfer functions between observed components 
of variation fields preserve amplitude, pnase and 
polarization information and have proved invaluable in 
the interpretation of linear magnetometer arrays where 
data from different arrays often have to be combined to 
define a geomagnetic induction anomaly (Schmucker, 1964 


and 1970a; Everett and Hyndman, 1967). 
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Following Schmucker (1970a), suppose the function 
Z(t) is linearly related to X(t), both normalized to zero 
mean value in the interval -T/2 < t < T/2 . Their relation 
shall be indenendent of time in this interval and thus 
expressible by a linear transfer function zy in the frequency 
domain, assuming T > ~ . If we denote the Fourier 
wants forms sO fer (t) aide Ut) ainseg. eaaheby COX )eand —€ CZ) 


respectively 
CUA Z y C(X) (2ee\5o) 


Furthermore 


Sis ; (2.16) 
SGN ELS. (gay = SE ew (2.17) 
pie es) (2.18) 
SG Oe 


* * 
Here C (X) and S (XZ) denote the comolex conjugates of C(X) 
and =s( 2%). S(%) 1s) the auto-power of X(t) and S(ZXxX) the 
cross-power between Z(t) and X(t). If Z(t) contains some 


uncorrelated part 6Z(t), its power is given by 
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since [zy apse is the power of the related part of 


The anomalous components of the total geomagnetic 
variation field are related to the normal components of 
the field and since Maxwell's equations are linear in 
the field variables, the anomalous parts are expressible 
as linear functions of the normal parts. These relations 


can be expressed in matrix form as 


C(x.) Xy Xy Xz C(x.) C(6X) 
Ce) ai Pac eat ye aay, crane | aera AC. (ctiay) PCS 2A 
Oy Zy Zy Zz C(Z_) Ci(Gs7Z,) 


with x, and x, and the like representing the anomalous 
and normal parts of the fields respectively. The 3x3 matrix 
is the conplexs transfer function matrix such that 

Zee Zy(u) + izy(v) (2.22) 
for example, with Zy(u) tne real part and zy(v) the 
imaginary part of Zy. LiewLaste Vector al eq... 2.cdcontveains 
the residuals (uncorrelated parts) of the field components. 
IneeMatti nwo, traister LuUnGtLions must bes determined ssuch chat 
the powers of the residuals S(6X), S(6Y) and S(6Z) are 


minimal. 
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The third scalar equation in the matrix 
equation 2.21 can be written as 


C(6Z) = C(Z.) - ZyC(X_) - Z re) Daz C(Z_) 


With the aid of this equation and eq. 2.19 the power 
Specenum Of GlOZ 1S obtained. ~On differentiation 


this yields 


aS\o2) o * 
azy(u) ~ 7 HOS Za) SG TN eae Bez Ol Key 
(20:24) 
CVa * = 
IZ y ies Tues Lyn G ~~ - C (5Z)+C(X, )} 


The condition of minimum power in the residuals 
implies that both derivatives in eq. 2.24 go to zero. 
This is achieved when the cross spectra between §Z(t) and 


the normal parts Ky? fe and Zh vanish, thus when 
S(SZX_) = S(5ZY,) = S(6ZZ_ ) = 0 (2m 5) 


be vthe cross spectra of CioZ) with C(x.) Ce) 


and C(Z, ) are formed using this condition and eq. 2.23 


we obtain the matrix equation 
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S(X,) S(Y Xx) S(Z_ Xx.) Zy S(Z,X,) 
SOX Y) St Size.) Zy = S(Z.Y,) (2772 6) 
S(X,Z.) Sine) S(Z_) Zo Sa) 


and the vertical transfer functions Zys Zys Zz divead 
particular frequency can be obtained by matrix inversion. 
Similarly the horizontal transfer functions Xys Xv» Xz 
and Yy> Yy> Yz can be obtained. Because of the practical 
problems involved in separating the normal and anomalous 
fields, a topic discussed in Section 1.2.4, this general 
form of transfer functions (eq. 2.26) has hardly been used 
in geomagnetic depth sounding. 

Reliable estimates of transfer functions require 
the use of a number of events of various polarization and 
if possible due to various current sources, which are 
combined to obtain average power and cross-power spectra. 
Schmucker (19790a) showed that when M events are used the 
residuals should not exceed ((M-5)/(M-3)}? as upper 
permissible limit in order to have meaningful coherences. 
This means power and cross spectra should be obtained from 
more-than “five events. 

In practice the major problem in computing transfer 
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discussed in the next cnapter was the first large two- 
dimensional array for which transfer functions were 
calcuhated. For that study and also for the subsequent 
study in 1972 the following aporoximations were made. 
Firstly, for each spectral term of each event an estimate 
of the normal field in X and Y was made by taking mean 
costne and sine Fourier coefficients over all stations. 
Furthermore, at the magnetic latitude of 
southern Africa the normal Z component in substorm fields 
is small. The magnetograms in Chapters 3 and 4 support 
the view that the Z variations are mainly anomalous. It 
is thus explicitly assumed that Zi = 0 and Z. = Z So 
that the terms in eq. 2.26 that involve Z vanish. 


Consequently we can solye eq. 2.26 explicitly for Zy and 


Zy to yield 
S(ZX_) oe) - S(ZY,) SO | 
AST CIS ene cl ae 
n n neon (2.27) 
SLT Shy) ees GE) ash Xess) 
ee = 
3) Sth meee Cle Me 


Parkinson (1959, 1962) and Wiese (1962) introduced 
the use of arrow representations to display the anomalous 
behaviour of the vertical magnetic field with position. 
Transfer functions can be used to form similar convenient 
induction vectors. With 7 and 5 UNL tT Vectors ine thes qeo- 


graphic north and east directions at a particular station, 
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the induction yectors at anv particular frequency are 


given as 


“oe > + 
ie = Zy(u)i + zy (u) J 
C252:00) 
oe ee t 
Wy *; Zy(v)i x Zy(v)4J 


The real parts of Zy and Zy give the in-phase 
responses or 2 to Ne and Ae and their imaginary parts 
the quadrature-phase responses. When plotted the in-phase 
vector points away from a conductor in which the current 
flows in phase with the normal field. This vector is 
usually reversed to be in the same sense as a Parkinson 
arrow so that it points toward a conductor in which tne 
current flows in phase with the normal field. The quadrature- 
phase vector points awav from a conductor tn which the 
induced current leads the normal field in phase. The 
direction of this vector is also sometimes reversed. 

In the Earth induced currents flow at a phase 
angle between 0 and 1/2 with the normal field and a given 
Conductor will contribute to both induction vectors, “As 
pointed out by Gough et al. (1974), in the extreme case of 
very high reactance, arising typically when the normal field 
does not penetrate through the conductor, the in-phase vector 
will dominate. At the other extreme, when the normal field 


penetrates easily, as through a thin surface conductive 
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layer, the inductive reactance will be low compared to 
the resistance and the quadrature-phase yector will 


dominate. 


2.6 Numerical modelling 


ANne-dimensional modelling of layered structures 
is discussed in some detail in Chapter 5 and will thus not 
be dealt with in this section. 

It is often possible to consider observed mag- 
netic variation anomalies as approximately two-dimensional 
on a local scale and also to consider the inducing field 
to be nearly uniform over the region of anomalous induction, 
especially in geomagnetic mid latitudes. Under such 
conditions the modelling of the observed variations by 
induction in two-dimensional conductivity structures is in 
most cases adequate. The most prominent exception is when 
the anomaly is of the current concentration type. However, 
when no complicating circumstances exist, the general 
interpretation of induction in the earth is then that of 
considering the local perturbations (in two dimensions) by 
POCadmealLeralecOnaulcuily it yims UnUC UU Gem Ojon bieCaeb dil Le TaN yO. 
induced currents flowing in the large scale stratified 
(one dimensional) structures remote from the anomalous 
region. 

For this situation let the Cartesian co-ordinate 
evstemeixsy5z2) be aligned such that the x-axis 1s) parallel 


to the strike of the anomaly and z > 0 as usual is 
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directed downwards. Maxwell's equations (eqs 1.1-1.7) 
can now be rewritten for the two dimensional case where 


all quantities are independent of x and (X,Y,Z) are the 


cartesian components of B. BY combining €qs 122 “and 1. 
and assuming as usual that displacement currents can be 
neglected one obtains 
on Ol 
oa eo Woe (2 
Oke 
bea = UO y We 
Ons 
oy: = 2 Wo ES (2 
By combining eqs 1.1 and 1.6 one obtains 
JE JE 
a0 = = = - 1wX G2 
EY 
S71: -1wY (es 
JE 
ee ary (ee 
dy 
where all variables have a time dependence of exp(iwt). 
The requirement that 2 = 0 has decoupled the 


0X 
xX components of E and B from the other components of 


Cache vector.esovtenat eqs 2.29 to 2.34 separate into, two 
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independent sietsaof equations. Equations 2.29, 2233 and 
2.34 involve only Es YeseZeWi We Regs 4o0 4.2.45 lecaldiezen 2 
involve only X, by and ES: BY combining theativst set 


and eliminating Y and Z one obtains 


, 

anes a Ey 

ee + ei Some’ (2.35) 
y FL 


where E is parallel to tne strike and is termed the 
E-POLARIZATION case. 


Elimination of EY and ES from the second set 


yields 
2 2 
eae ae ‘Werte (2.36) 
MY Ze 


which is termed the H-POLARIZATION case. 

The problem of adapting these diffusion equations 
LOmdaMunerical solution of observed anomaliles nas been 
extensively studied in recent years (e.g. Weaver, 1963; 
Wright, 1969; Hadden and Swift, 1969: Jones and Price, 
iO Ole Oe) ee Wie been Oy le ONS mama ras COCO) Me Coggon:, 
1971: Pascoe and Jones, 1972; Reddy and Rankin, 1973). 


Equations 2.35 and 2°36 can be written as 
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where ae toeeither by or Ke Nim = Ome and the subscript m 


identifies the different units within the model. In 
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geomagnetic depth sounding only the E-polarization case 
is of relevance because the H-polarization case does not 
produce anomalies at the Earth's surface. The relevant 
boundary conditions have been the subject of some 
discussion ‘(ct. Weaver, 1963 Vo Jones and Price. 1970; 
Hermance, 1972 v. Price and Jones, 1972). The computer 
formulation used in this thesis is based on the algorithm 
of Jones and Price (1970) as programmed by Jones and 
Pascoe (1971) and Pascoe and Jones (1972). The boundary 
conditions as qiven by Pascoe and Jones are used in the 
computer program. 


The Jones and Price (1970) method is a finite 


difference metiod and with i = i a ig. po PO Gne ee oy 
separates into 
2 @ é 
Vet 7 Um Om 
(253:3)) 
2 rye 
Yo 9m = Un tm 


These two equations are replaced by their finite 
difference equivalents which are then solved at a rectangular 
mesh of grid points under the condition that at the 
boundaries between mesties, the tangential magnetic and 
electric field components are continuous. By the use of 
the Gauss-Seidel iterative method, the solutions at the 


mesh points are iterated until stable convergence is obtained. 
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A linearly polarized external field with 
uniform horizontal component at a large distance from 
the main conductor is applied to the system as a source 
field, and the total field is set to zero at the greatest 


distance into the main conductor reached by the grid. 
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CHAPTER ltd 
THE 1971 ARRAY STUDY IN SOUTH AFRICA 


Se eelntroduction 


A triangular array of magnetometers was operated 
in central South Africa from August 3] to October 27. 
1971 by the Geophysics Division of the National Physical 
Research wabhOvdatonry 10 @EnemG son kane Cretormias. South 
Africa. The 24 three-component magnetometers were of the 
type described by Gough and Reitzel (1967) and were on 
loan from the University of Alberta, Edmonton, Canada. 

Rie TOCatvone Ot  thesarray 1s SHOWN in) Fide. 3.4 
Inarelavion wo the tectonivce provinces: and to the conti- 
nental edge. The study was designed largely to look for 
gepossible contrast tm electrical conductivity of the 
crust and upper mantle between the Kaapvaal craton with a 
typical basement age of 2600 Myr and the 1000 Myr old 
Namaqualand-Natal mobile belt (Nicolaysen and Buraer, 
1965). The distribution and main features of these tec- 
Lover pDRovinces were discussed 1 section 1.3.) In section 
1.4 it was pointed out that the resistive portion of the 
crust in the Kaapvaal craton as determined by ultra-deep 
Schlumberger soundings seems to be in general more conduc- 
tive than the same unit in the Namaqualand-Natal mobile 


Hebel tewas: realized that the crustal, conductivity con- 
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trast between the two regions, which differ by a factor 
Of about Web. Was) not enough to cause a detectable induc- 
tion anomaly, but at that stage we had no knowledge about 
possible conductivity Hilt Peences in the upper mantle. 
The boundary between the Kaapvaal craton and the 
Namaqualand-Natal belt is indicated in Fig. 3.1 and the 
arcay screacdted at. 

The field procedure similar to that described 
excellently by Camfield (1973) will not be repeated here. 
The instruments were serviced at about 21 day intervals. 

Station code names are given in Fig. 3.1 and 
the geographic co-ordinates and magnetic declinations 
at the Stations are given in fable 3.1. [he work discussed 
(Hetiismendpteratace Deen tne sSUbIeCtOhed paper by «Gougn, 


DeBeer and Van Zijl (1973). 


3.2 Magnetograms 


Three periods of geomagnetic disturbance were 
selected for analysis. These periods were August 31, 
18.30-22.30 G.M.T., September 17, 20.00-September 18, 
06.00 G.M.T. and September 26, 15.30-September 27, 03.30. 
Mneerecorainguerticlencies 10n the three vents were 
respectuve lyec0u,a7 47 and 792. che main causes of loss 
Orerecord lay in poor operation of the Stepping motor 


camera drives and in traces going off scale or becoming 
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Table 3.1] 


Coordinates and magnetic declination of stations 
in the 1971 array 


Latitude Longitude 
Degrees Degrees Magnetic 
Station and min and min declination 
Name Code South Bast (degrees) 

Prince Albert PAR Sve Sey" ale On Sees 
Road 

Fraserburg FRA 31°40" Cann Oo -21.9 
Beaufort West BEA 32. 09% (aa aoe -21.9 
Biesiespoort ie Sled oy: Coal en =2ie/ 
Mynfontein MYN S000" ESAS a -21.1 
Petrusburg Beat SO a= CA 35 -20.7 
Koffiefontein KOF EAS 2A Sie -20.1 
Dealesville DEA ZO Oe ASPEN) -19.8 
SCuMI aes dri ht Sol 23, 42 24°04! -19./7 
Hiag's Hope HIG EMS oe Cea e =20 63 
Prieska PRI a as Pie Oe -20.9 
Vanwyksvlei VWY 3.02.38 ANS kes = (aa ee 8 
Tontelbos TON SUC by ZO eal meyers) 
Granaatboskolk GRA 30°02" Tobe =.) 
Pofadder POF 29°04! 19°34! -20.4 
Kakamas KKM ZG. 4 Om 203 oe -19.5 
Ventersdorp VEN 2 0e Gn 26 54» -17.8 
Christiana CHR AU NS 25° 0G! -18.9 
Vryburg VRY ci as ZOU: -18.5 
Kuruman KUR Cieeae Cat le -19.2 
011i fantshoek Oia fase OTE (Bie Taine =O 95 
Vredefort URE No data were obtained from these 
Upington UPI stations due to instrument problems 
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too faint jon the film. The selected data sections were 
digitized and reduced for analysis as described in Section 
2.1. The sets of normalized time series representing the 
three orthogonal components in the geographic coordinate 
system (X northward, Y eastward, Z downward positive) 
were plotted in stacks by means of a calcomp plotter. 
These magnetograms are illustrated in Fig. 3.2 for the 
SUDSTORMPEVENG OTeAUGUSE SI, 91h =hiders.3 forethe- distur 
Dancesevenu1 foeuLemberm 1 / sande loo andain st! ioe ss.4) tor 
the event of September 26 and 27. The three ‘lines' of 
stations run from the three stations nearest the north- 
east corner of the array south-westward, Line 1 being the 
north-west side and Line 3 the south-east side of the 
triangle (Fig. 3.1). In each stack of magnetograms the 
sequence from top to bottom runs from north-east to 
south-west along the line. The time marks are in 
Universal time (UT) and the local midnight in the centre 
of the array is around 2230 UT. Scale bars relate the 
amplitudes to gammas (1 y = VOee et enim 07 oersted). 
Note the different scales for the vertical and horizontal 
components. 

Visual inspection of the sets of magnetograms 
in Figs 3.2, 3.3 and 3.4 indicates that the three events 
have the following common features due to local anomalous 


variation fields superimposed on the normal part. The 
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Figure 3.2 Magnetograms of a substorm 18.30-22.30 G.M.T. on 
AUGUS tEO le toy Im ln ea chen ine Sthe es tations acun 
north-east to south-west down the stack. 

Components are X northward, Y eastward and Z 
downward. The unit is the nanotesla (gamma). 
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Figure 3.3 Magnetograms of a magnetic disturbance event 
September 17, 20.00-September 18, 06.00 G.M.T. 
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Figure 3.4a Magnetograms of a magnetic disturbance event 
September 26, 15.30-September 27, 03.30 G.M oie 
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Figure 3.4b Magnetograms of a magnetic disturbance event 
September 26, 15.30-Sentember 27, 03.30 G.M.T. 
Line 2. 
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Fiqure 3.4c Maqnetograms of a magnetic disturbance event 
September 26, 15.30-September 27, 03.30 G.M.T., 
Leno. 
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magnetograms for each component are rather constant along 
Dine tit, CXGept ror a slight increase im «the amplitude of 
Z to the south-west. In Line 2 there is a marked increase 
in Z and a smaller increase in X south-westward. Line 3 
has the same trend, except that the Z component shows a 
large enhancement at Beaufort West (BEA) and is greatly 
reduced at Prince Albert Road (PAR). In Fig. 3.2 the 
P traces at tnese two stations show this difference inva 
Striking way for periods of about 30 minutes and 2 hours. 
The X component increases steadily to PAR at the south- 
western end of Line 3. Comparison of traces along the 
FOWS@ IN Soc. oo ONG. sede SHOWS —that Land \ i1ncrease 
to the south-east, and it is clear that the main effect is 
a southward increase of X and, except at PAR, of Z. 

The waveform Z(t) bears a strona resemblance 
to X(t) ab the more southerly stations, and still resembles 
X(t), passed through a high-cut filter, even at the 
north-east of the array. There is no resemblance of Z(t) 
to Y(t). The amplitude variation of Z and the resemblance 
of its waveform to X suggests immediately that much of the 
Z variation Pea is an anomalous field resulting from 
induction by X in a conductor near PAR at the south-west 
COnneL Ion che drrdy.  Ihis conductor also. produces <a marked 
attenuation of high-frequency variations in Z at PAR 


reminiscent of the attenuation of Z west of the Northern 
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Rocky MounitaincavCamTield et ale. 1971). Thesattenuat on 


SUGGESTS tna teen tse above the conductor. 


Seo apse On TroOuUrter spectral components 


Sequences were selected for Fourier transforma- 
tion from each of the three data sets described in the 
previous section. These data sets were August 31, 18.30- 
eeo0 eG. Mela. September 1/5 22.00=September 18. 01650 G.MoT. 
SED CEMDEN S2Og hGwo0-27,50 GeNe tl and September 26, lo. 30= 
September 27.03.30 0G.N, i. In each case -anblatude spectra 
lvkesthosesshown ins Rugs 325, 3.6 "and 3.7 “for «the 
four-hour long sequences were plotted. The periods used 
in constructing the contour maps were chosen near spectral 
pedks 10y and Zato gev the optimum signal to nore: ratao7 
Contoured maps of amplitudes and phases at the selected 
periods were prepared; in all ten sets each containing 
Six maps were drawn. The techniques are described by 
Reitzel et al. (1970) and Camfield et al. (1971). The 
period range covered is 24 minutes to 293 minutes. The 
12 hour long data sets for September 26-27 were transformed 
toe nrovide then mans: at this last period.  Pherpolari zation 
of the horizontal field was determined using the theory 
descmiped Minuscection) 2.4.5 (he polartzabion ellipses were 
always elongated and often nearly linear, with major axes 


rangingerrom N5o-W) to N2@°E for the ten periods mapped. 
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Fiqure 3.5 Fourier amplitude spectra in the period range 
20-240" min’ for 10030-22030 0G. Me eeeerG Uc teas) 
for the stations VEN @andesBEens 
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Figure 3.6 Fourier amplitude spectra in the period range 
20-240 min for September 17, 22.00-September 18, 
Ole SOSG. tel. TOY the Sta crons .VRY and: BEA: 
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Figure 3.7 Fourier amplitude spectra in the period range 
20=240 min for September, 26, 16230=22.30) GoMiT. 
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Typical sets of Fourier spectral component maps 
are shown Ih Figs 358,°3.9° and 3/10; “These maps cover 
the period range 24-120 minutes and polarization azimuth 
range N21°W to NI17°E. All of them show a large anomaly 
in Z with a maximum near BEA and lower values at PAR, as 
already noted from the magnetograms. All maps also exhibit 
a southward increase of X suagestive of a maximum in this 
component a short distance south of PAR. The amplitude 
maps of X and Z are thus qualitatively consistent with 
the hypothesis that X induces current in an approximately 
east-west striking conductor just south of PAR and that 
Z is mainly an anomalous field produced by this current. 
This hypothesis is strongly supported by the phase 
relationships. In all maps the phase of Z is very close 
tO-tiate Of Keatestatrvons suchas BEA Boe and BRA whitch 
are near the amplitude maximum in Z. This phase relation 
ise shown Ine rigs §o.8=3. 102 ~In fable Gx2 the phase 
difference (Z-X) is shown at BEA, at the maximum of the 
ZmAnOMalyeeanduateOl lata: tOmthesmnOrcn. me Lumtseannteres tind 
to note that the largest phase differences between Z and 
YROCCCUL MOT atiemil is tetwo pe niOdS -dtanire COP 10 f sCiter cabies 
these have. the largest content of Y in the horizontal tield 
and north-west polarization. 

A small maximum in X lies near the stations 


OLI, PRI and PET in Fig. 3.8 with an associated maximum in 
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X AMPLITUDE X PHASE 


Y AMPLITUDE ; Y PHASE 


Z PHASE 


Fiqure 3.8 Fourier transform amplitudes (nT/c/min x 10) and 
ohases (min) at period 24 min for 18.39-22. 30 
Gollele es RUdUS LD 31. “ES timates for tes tandeard 
deviations in the X, Y and Z components are 
respectively 18, 17 and 14 nT/c/min for the 
amoliivaes and <5. 2 atcel-5 amin. for che 
phases. The polarization ellipse for the 
horizontal field is shown. 
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X AMPLITUDE X PHASE 


Y AMPLITUDE Y PHASE 


Eragure 3.9. Fourier transform amplitudes (nbi/c/min =< 10) and 
phases (min) at period 51 min» for 18. 30-22. 30 
G-M.7., September 26. Estimates for the standard 
deviations in the X, Y and Z components are 
respectively. 20)., LoxandalS ni/c/min.for the 
amplitudes and i, | and [lso=min for uthesphases: 
The essentially linear polarization of the 
horizontal +1410 1S. indicated. 
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Figure 3.10 Fourier transform amplitudes (nT/c/min x 10°} 
and phases (min) at period 120 min for 18.30- 
Cems eGen lee AUGUS tC Slee estimates «hor stime 
Standard deviations in the X, Y and Z components 
are respectively 43, 23 and 25 ni/jc/min for the 
amplitudes and 1.5, 1 and 2 min for the phases. 
The polarization ellipse for the horizontal 
field is shown. 
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Le JUSTONOriIne= On tnis. Similar anomalies are found in 
another map (not shown) at period 43 minutes for August 31, 
rorownich the horizontal field is also polarized NW-SE. 
The X anomaly coincides with the boundary between the 
Kaapvaal craton and the Namaqualand-Natal mobile belt. 
Since the current system giving rise to the 
large, southern anomaly is beyond the corner of the array, 
contour maps are of little use in determining the geometry 
of the current system. Even the question whether the 
CONCUCTONTISe Straight orscurved is open. . iranster 
functions were therefore computed to try to answer this 


question. 


524 Transfer functions 


The centre of this array was situated at magnetic 
latitude 40°S. The magnetic latitude (uy) is derived from 
theamaqnetac inclination) (1) according: to <2 svanu.— 9 Gane le 
hte lonaicudeme ck whi cherUns., through tie centre of ithe 
array the geomagnetic and magnetic equators differ by about 
Otel he mogneticmlatitude si ssquotedubecause the miiat 1 iu 
dinal variation in geomagnetic activity is determined by 
the magnetic and not the geomagnetic equator (Rastogi, )19o2). 

At the magnetic latitudes where this array was 
Situated, the normal Z component in substorm fields as 


Sialioe the maqnetourams Of Figs 322-354, the sphase 
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relations ine Table 3.2 and the amplitude spectra for the 
Sat ONSeM a theamnortherne part of the array (rigs §3e5=3.7.) 
Support the view that the Z variations are mainly 
anomalous, and related to induction by X in a conductor 
SOUThDOn the larrayeseVenmats thernorthern stetions. Lt. is 
therefore reasonable to calculate transfer functions from 
components Ls and ‘ OT whe Nornal field to Z, where 14 
is assumed that ZL. = 92. Anes timateno; ther normal =fieid 
in X and Y was made by taking for each spectral term of 
each event mean cosine and sine Fourier coefficients over 
alties tat Ons. She transter funetians Zy and Zy were 
Caiculated using the relations given in eq. 2.27. 

(Nemo anSie nun CE TOMSs 7, .aaiids:Z 


X if 
quantities whose real parts give the in-phase responses 


are complex 


OuhenZe GLO X and ee while their imaginary parts give the 
quadrature-phase responses. The real parts of zy and zy 
at a specific station were combined to yield an in-phase 
induction vector which was as usual reversed to be in the 
Sames sense as a Parkinson arrow (Parkinson, 1959, 1962). 
The vector thus points toward a good conductor in which 
Clrrente tlowseanephase With the normal field. Une 
Timdidanarys Parcs of zy and ay were similarly combined to 
give a quadrature-phase induction vector which points away 
from a conductor in which current leads the normal field 


in phase. 
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inmoection c.55. on the scaleulation of transter 
TUNGUIOUS  WeuWwas -pOINted. out that each cross= and 
autopower must be a mean value from a number of variation 
events with well distributed polarizations. Seven short 
events each of 1290 minutes duration were therefore 
selected from the three disturbance sequences. 

The in-phase and quadrature-phase induction 
vectors repnresenting zy and zy derived from seven short 
events are shown in Fig. 3.11 at the fourteen stations 
which recorded all seven short events in three components. 
Four periods are represented. The in-phase vectors indicate 
the presence of a good conductor south-southwest of the 
array, and the parallel orientation of these vectors 
indicates an approximately straight current confiauration. 
The large difference between the in-phase vectors at PAR 
and at the stations BEA and FRA indicates that PAR is 
close to, the conductor. 

In Section 3.3 it was remarked that the attenua- 
tion of short-period Z in the magnetograms at PAR suggested 
tietethis station was’ apove the conductor, and the Steep 
rise in X amplitude as PAR is approached from the north 
imcdce bese that the current CONTriguration 1s) just south of 
PAR wIihese facts Strongly stlgdest that a large conductive 
pOdveltesmunder the cape Folded Belli (Fig. 3.1) The 


continental edge is too far away to produce the steep anomalies 
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In-phase (solid) and quadrature Puget (broken) 
induction vectors representing the transfer 


functions zy and Zz) at four periods, f-omes cua 
2-hour variation events. 
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in Zand X-indicated in. Fids 3.8-3.11. “Although, 

the in-phase induction vectors are not perpendicular to 
the continental edge (Figs 3.1 and 3.11) they are parallel 
Toy tresRarkinsonsdnrOW ofor Hermaniisa Grids sal). a 
magnetic observatory on the coast some 100 km west of 
ThessouLtnern tipeot eAirica (Parkinson, 1962), 

An important point is the continued response of 
the conductive body as the period of the fields increases 
to 64 and 128 minutes. This frequency response is 
consistent with a deep conductor of large self-inductance, 
and so with an upper mantle structure. The in-phase 
KesponsSe. On =/. (tou Indicated aby -Fldmessali) and Fab lence 
are consistent with a conductor in the upper mantle. 

The southern corner of the array lies across 
the deep Karroo sedimentary basin. The effect of this 
Wivueeven GiscuSssedo a Nese Clon: cc Orn = NmetGie conaltlwletecan 
be seen that the quadrature-phase induction vectors at 
BEA, FRA and PAR show the effects of the sedimentary basin. 
At these stations the quadrature-phase vectors are larae 
at short periods and are reduced to magnitudes like 
those elsewhere in the array at periods 64 and 128 
minutes. At the shorter periods the quadrature-phase 
vectors tend to point away from the deen Karroo basin. 
No doubt the in-phase induction vectors at PAR are 
shortened, at the shorter periods, by the effect of the 


sediments to the north of this station. 
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3.5 Normalized anomalous fields 


Anomalous fields cover most of the array. Northern 
Stations, as far as possible from the southern anomaly, 
were arbitrarily desiaqnated as representing normal field 
components. For X and Y the mean amplitude from KUR, CHR 
and VRY was used to define the assumed inducing field over 
the array and, as residuals, the anomalous horizontal 
components Mee i at other stations. For Z the mean 
amplitude from VRY, CHR and VEN was deducted from Z at 
each other station to give an estimate of the anomalous 
vertical field Zoe A profile was chosen approximately 
parallel to the in-phase induction vectors: it is shown 
Tega le mE rOlmcOn FLOUR EmMab Sess Mit lavrestOu i 1OSen Gao lO 
and the normal field estimates just specified, normalized 
anomalous fields ate and He te were estimated along the 
profile. Here H is the horizontal component along the 
profile. Normalized anomalous fields for six periods 
from three events are shown in Figs 3.12 and 3.13. Four 
other anomaly profiles, not shown, resemble those of 
Fie ole which Shows sthesresponse: tO Tlelds polarized 
nearly north-south. The more complicated and larger 
anomalies of Fig. 3.13 arise when the normal field 
polarization is close to north-west-south-east. 

The normalized anomalous vertical fields i a 


associated with north-south polarized normal fields (Fia. 
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Figure 3.12 Normalized anomalous variation field components 
Zi/H. and ue al along the section line AB, 


Fig. 3.1, for periods 24 min (event of September 
17, 1971), 51 min, 102 min and 293 min (evete 

of September 26, 1971). The polarization 
indicated in Fig. 3.9 is representative (east 
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Fiqure 3.13 Normalized anomalous components Bsa ale and Hoe 


along AB (Fig. 3.1) at periods 24 min and 43 mam 
in the event of August 31, 1971. The polarizatime 
in Fiq.. 3.8 is representative (west of norte 
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3.12) show maxima which decrease as the period increases. 
At PAR ay TS smalivat 1 = 24 mineand riseseas. 7 
mcereases, suqqesting a southward shift) of the current 
system with increasing depth. The Karroo sedimentary 
Dasine Maye produce sone or all eorethasSeetrect.) sat short 
periods ill rises steeply towards the southern corner 
Ons tNe array, Duty this anomaly flattens outrat  oeriods 
over 100 minutes. Qualitatively the horizontal component 
anomaly also suggests a southward shift of the current 
System with increase of T and of the depth of the induced 
currents. 

When the normal field is polarized north-west- 
SOUth=east (Fig. 3213). both og and Hol develop 
extra minima and become much laraer than the corresponding 
anomalies with north-south normal fields. This feature 
of the response is not understood on the present limited 
information from an array which does not include the 
anomalous currents. 

Both» anomalies in Figs. 3-13. and those at | = 24 
ime nde) =e olemi ne iki. eseul2eare TOO mlardes tos pe 
compatible with induction in a conductor of two- 
dimensional geometry. The very large amplitudes must 
be associated with a three-dimensional configuration of 
COMGUCTON SINC wINGUuCt1 On OCrOWS On LG .eoe Hin end so 


Support to the presence of a curved conductor. One 
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POSSsibility- 1s that the anomaly is of the “current concen- 
tration’ type, with induction in a large region feeding 
a highty-conducting body. A second possibility is that 
constructive interference is present between the fields of 
two Or more current systems. One could be in the Karroo 
basin and the other under the Cape Folded Belt. 

Whatever the inductive mechanism, a maximum 
depth can usually be estimated from the half-width of 
an elongated anomaly defined as the width of the Ha 
anomaly at half height. For a line current the half- 
width is twice the depth of the current, and the Z 
extrema are a half-width apart. Unfortunately the 
DGOnl es 01 | lGS= esuleaancd culo doOenoL extend far enough 
to the south to allow the half-width to be estimated. 

Normalized anomalies at periods 24 minutes and 
102 minutes are repeated in Fig. 3.14 with a section 
showing the basement topography and the resistivity 
structure in the Karroo sedimentary basin. Average 
resistivities are shown as deduced from various deep 
structural electrical investigations whereas thickness 
information was obtained from deep drillings and seismic 
Yeflection work (Haughton, 1969: Fatti and Du doit, 1970 


and Winter and Venter, 1970). 
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3.6 Two-dimensional numerical modelling 


Despite the fact that the arvay covered only 
about half the width of the anomaly, two-dimensional 
modelling was attempted to attain a better understanding 
of the anomaly. The two-dimensional model calculations 
have been made by means of the algorithm of Jones and 
Price (1970) as programmed by Jones and Pascoe (1971) 
and Pascoe and Jones (1972). This method has been 
desemiped in section 2.6. The models all include the 
Neawesuinace: Vesustrivity date Shown. (i Fides ol4. 
Thirty-five different models were calculated. It was 
first established that the sedimentary basin alone 
could not be responsible for the anomaly, even when 
conductivities were scaled up Dy a factor of ten to 
allow for a possible current concentration effect. The 
computed anomalies were too far to the north and too 
narrow in relation to the observed profiles Zfu, and 
Ho/H,- Various conductors of rectangular and more 
complicated sections were added to the sedimentary model, 
at e@rtierwcrustal OM Upper Mantle depths.  sResults. for a 
crustal and a subcrustal conductive structure are 
respective|yusnown in Figs j32l5 and 3716.) The resis- 


tivity data for the “normal" crust away from the anomalous 


bodies correspond to that given by Van Zijl et ale at-957 C3) 
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Results of two-dimensional model calculations 
for a compli cated “crustal structure: with 
resistivity of 19 ohm.m. The normalized 
anomalous components (ESA and nisi Pists along 


ABs Phasstile ate pertod of 24° mitiiaremfor the 
event of September 17, 1971 and at period 120 
Mite kOe ecie event eon “NUGUS ts 3. 1971 gee hie 
resistivities in the Karroo sediments are 
CHOSercindicatedeih waidss o6 14 


> f 
4 
on 1 
‘c 
™" 
i 
: a 
’ : 
’ 


- 
: . - 
‘ Av 
~~ s' 
hy 
7 : ve toogme . 
——- os ‘ fe e 
\ ’ 
\ i — FY 
/ 
/ 
J ~~ = 
a ty 
! 
a Pr ' 
1 \ 
a4, .j@ : 
| i 
sitiia a ~=y, 
7 ' 
1 
af) 
aoe 
_ _ 
- 


pooh wa le? 


19 ty 


tg 


sqitinvis 158 neta | 
ete gt cre 


PY ai 
ips 


bots 


gt pang | 


- 
L) 
! 
4 
« @ 
¥ ar « 
qe ee 


ve 


7 ww ae 
labo fenthe dan Hy ow 4 6 


1 aor eee 


2g 
Hy 


Braure: 3.16 


+ 


OBSERVED ° 
COMPUTED -=——- 


Results of two-dimensional model calculations 
fOr a Sinple suberustalestructure: with 

Lesis MiVUty. Od coomOlMmeinle) lie Zo/H, and Tye 
data are the same as those shown in Fig. 3.15. 
The resistivities in the Karroo sediments were 
decreased to 1/10 of tne values shown in 
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for the Kaapvaal craton and Namaqualand-Natal belt. These 
Cota diiten slightly from those given by Van Zig] and 
Joubert (1975) as discussed in Section 1.4, but were 
the only available in 1972 when these calculations were 
made. Since the differences only involve the very resis- 
tive strata, tney will have little or no effect on the 
COSC Ss 

It proved possible to model either Zi/t or 
sey ane but not both, all models making CHE a) Bees) too 
large in relation to the observations. Tits pron Vem as 
Vevyewelm Illustrated in t1gs. 3,15. end 3.10.) From 
these figures it is clear that the frequency response 
also could not be matched by the model calculations. 
This was true for all models computed. In this respect 
Porath et al. (1971) encountered similar difficulty in 
trying to model the North American Central Plains anomaly, 
which is known to involve the current concentration effect. 
It was found that most models based on anomalous bodies 
TP neacrus tp rocdluce ae and ean anomalies that are 
too narrow in relation to the observed profiles. Anomalous 
conductive structures in the mantle gave better results 
ime this respect. The best restiits were obtained Using 
models incorporating enhanced conductivities in the Karroo 
and Cape sediments, as well as mantle structures. As was 


remarked in the last section, the present anomaly is too 
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large to be produced by induction in a two-dimensional 
Structure, so that two-dimensional model calculations can 
avebest locate the structure on some assumption. as to the 
amplifying effect of current concentration, constructive 


interference or some other cause. 


Or) OS CUS SALON 


Although the attempts to model the magnetic 
induction anomaly were unsuccessful, it is clear that a 
hag Ore Structure, sInvolving conductivity contrasts: of 
at least two orders of magnitude, underlies the Cape Fold 
belt and the deep part of the Karroo Basin. The Geophysics 
WeVelsaLOnN On eLNeu Cush. alee CeLOllewe SOU UNM Tr lal olsand © 
present constructing an array of 26 Gough-Reitzel magneto- 
meters. This array will be used to extend coverage towards 
the coast. Quantitative interpretation will be attempted 
when the anomaly is more completely known. The qeoloay 
and other geophysical parameters in the array area will 
be discussed in Chapter 6. 

The small anomaly near the stations OLI, PRI and 
PET for horizontal fields polarized NW-SE could be 
associated with the boundary between the craton and the 
mobile belt. This boundary seems to coincide with the 
Doornberg fault zone, a fault in the granite-gneiss 


basement which seems to be rejuvenated in the overlying 
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Karroo sediments (Olivier, 1972). The anomaly cannot, 
however be due to local induction, because it only 
shows up for the horizontal field polarized almost 
parallel to the fault zone. This anomaly cannot be 
explained on the basis of this array study alone, but 
will be discussed in Chapter 6 in the light of information 
oOnetie los cosarray study  Lunthere moOrtunr 

No indication has been found of any conduc- 
tivity contrast between the Kaapvaal craton and the 
Namaqualand-Natal mobile belt. Any compositional effect 
On cle cOndUCcCEIVItTyY OF the crust or mantle in these 
PEGLONLTG UNITS 1S) toO smalie to be detected by the present 


Study. 
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CHAPTER IV 


THE 1972 ARRAY STUDY IN SOUTH-WEST AFRICA 
BOTSWANA AND RHODESIA 


aol lnc roduetton 


An array of 25 three-comvonent maqnetometers 
operated in South-West Africa, Botswana and north-western 
Rnodesta from December 13, 1971 to February 2, 1972. The 
POcetnOne Ot Le sartay 1S Siow In Fics 4. te. All thie 
instruments but one were of a type described by Gough and 
RPeitzel (1967). The instrument near Tsumeb (TSU) is a 
permanent recording station of tne Magnetic Observatory 
of the C.S.I.R. and is of a type described by Scheepers 
(HOT) Table “4.1 gives the Tames and coordinates of 
all stations in the array as well as the magnetic declina- 
POM ea cee aces tation. 

HMesDOSiELOne OF LNG saraeay 1S) SilOW ied) Ghnlc) meds 
in relation to the tectonic features and some earthquake 
enicentres in the area. Unconsolidated Kalahari sediments 
Vamdinid fades ron dentiary CO Recentacovar Most of 
BGtsewand and the eastern part of South-West Atrica so 
that very little is known about the underlying solid rock 
geology and tectonic features. Faults of Karroo to 
Pecent age in South-West Africa and Rhodesia are shown 
(Korn and Martin, 1951 and Provisional Geological Map of 


Rhodesta, 1971). 
leh 
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Coordinates and magnetic declination of stations 
in 1972 array 


Latitude Longitude 

Degree Degree Magnetic 

and min and min Declination 
Station Name Code South South (Degrees ) 
Eiseb EES aly Sie! 19°46! ~ahera 
Ntjituo OTJ OA 28 hoses. -15.4 
Rundu RUN livia o 5s 19°45! = (3x0 
Namutoni NAM Lena 0 16°54" -15.4 
Okaukuejo OKA 19°00" eye -16.] 
Gwe ta GWE 2034s fay ANion = 3m0 
Aranos ARA 2355 20°00" =l7 96 
Makalamabedi MAK ZO lls 2 Shee a =1.3 20 
Outjo QUT ZO lee 16o20) 0% -16.8 
Gobabis GOB VAT ae 1924" -16.9 
Happy Valley HAP ag RA ONO 20.0" -15.9 
Lone Tree LON 22°45" col 00: = 6.473 
Sukses SUK (a ae ae spree -16.9 
Sehitwa SEH CO cs ie Canin -14.2 
Nokaneng NOK oss eA -13.8 
Mariental MAR 24°30' S56" -18.9 
Rehoboth REH 23.20% W038. -18.7 
Dett DET LS es7~ AO Se Sees 
Windhoek WIN (OSATE IZ 206s -17.9 
Sepopa SRP. 18°44! AAA WS pe = 3a 
Mukwe MUK 18°04! Calh cale -13.3 
Kwai River KWA ES 0.9. Zoe as -13.0 
Lodqe 
Kavimba KAV ken Oras 24°36: -12.0 
WiC LOT aeta'l. Ss VEG BST ies Vey Shore -11.5 
Tsumeb TSU 1202. leOne 3 OF: -15.4 
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The major linear features in Botswana as outlined by 
Reeves (1972) are presented. Epicentres for earth 
tremors catalogued by Fernandez (1972) for 1971, 
Fernandez and Guzman (1973) for 1972 and Rhodesia 
Meteorological Services (1973) for 1973 are indicated. 
Paeves discussed the Botswana seismicity and concluded 
that the epicentres form two distinct snatial populations - 
that of the Okavango Delta and that of the Central 
Kalahari. These two seismically active areas Reeves 
believes to be separated by the Ghanzi Ridge, a broad 
upwarp of little explored Palaeozoic sediments isoclinally 
folded along a north-east trending axis. Outcrops of 
rocks in the Damara Geosyncline are also shown. 

The array was located to look for possible 
Contrasts in electrical conductivity, or the crust ana 
upper mantle associated with the Damara geosyncline 
and the seismicity in Botswana. A deep direct current 
electrical sounding (Schlumberger array) that was centred 
about halfway between SUK and WIN (Fig. 4.1) indicated 
“ido the resistive part of the crust. in the Damara 
Geosyncline has a transverse resistance (thickness - 
resistivity product) an order of magnitude less than that 
of the crust in the Kaapvaal Craton and Namaqualand 
Metamorphic Complex (VaneZi gi, Higomand=ae Bel locas 31970. 


Woneeziil and Joubert, 1975). 
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Since the magnetometers require film changes 
at three- weekly intervals the siting of the instruments 
was largely dictated by the few negotiable roads in 


the study area. 


4.2 Magnetograms 


From the geomagnetic data recorded over the 
nine weeks of operation, three periods of geomagnetic 
disturbance were selected for analysis. These periods 
were January 16, 14.00-23.00 G.M.7., January 21; 11.30- 
anulany 2258 04 so0m Geol. and. wantany aco. Wy 4esUece. sOe GeMihs 
Magnetograms for these events are presented in Figs 4.3- 
4.5. The H and D components (magnetic coordinates) were 
transformed to X and Y components (geographic coordinates) 
because of a 7.5° change in declination across the array. 
Stack let rids “as-455 portrays. the, data of tne 
westernmost stations from OKA in the north to MAR in 
the South stack 2 the data of the adjacent) “lines oF 
magnetometers from NAM in the north to ARA in the south, 
Stack 3 the data of the next ‘line’ from RUN to LON and 
Stack 4 the data of the group of six easternmost magneto- 
Mecvers:. 

The most prominent feature in these magnetograms 
7s the reversal in the vertical variation field between 


OUT and SUK in Stack 1. The reversal in Z occurs in less 
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Figure 4.3 Magnetograms of a substorm 14.00-23.00 G.M.T., 
January 16, 1972. The stations are represented 
from north to south? in) each stack as combonenres 
are X northward, Y eastward and Z downward. 
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Figure 4.4a Magnetograms of an event January 21, 11.30- 
Januaryecc, U4. G00Gs Melemencoib a Grose le 
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Figure 4.4b Magnetograms of an event January 21, 11.30- 
danuany #4c. 04330 9Ge Melee track ace 
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Figure 4.4c Maanetoarams of an event January 21, 11.39- 
Januaryic22., 204 730 @ GEM ier toc Kaos 
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Figure 4.5 Maqnetograms of a magnetic disturbance 14.30- 
22.30 GoM.l 2, January 2onw Losec seuiesme Ven 
was recorded with more than 90% efficiency. 
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than 100 km. The waveform Z(t) at OUT bears a strong 
resemblance to X(t) while the short period variations at 
SUKeseem to bey 180°. out of phase with those of OUT. The 
waveform Z(t) bears a strong resemblance to X(t) at the 
more western stations, and still resembles X(t) passed 
through a highcut filter even at the east of the array. 
IMS iss perricivarly welear in Figs 4.3.0 There 71s no 
resemblance of Z(t) to Y(t). The reversal in the 
vertical variation field in Z and the resemblance of 

its waveform to X suggest immediately that much of the 

Z variation field is an anomalous field resulting from 
induction in an east-west trending conductor between 

OUT and SUK: 

Other features in the magnetograms are the 
enhancement of the X variation field at SUK and EIS 
indicating that the maximum concentration of induced 
Current must be nearly beneath these stations, and the 
small amplitudes of the normal vertical force field in 
Botswana. Furthermore the high frequency variations in 
Z at EIS, LON, SEH and KWA show a marked attenuation 
which suggests that these stations are above conductive 
SurucLUres:. 

A reversal in Z is also observed between TSU 
and GOB, between NOK and HAP and between SEP and GWE. 


It is clear that reversals for short periods become phase 
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siimtseot Vess than 190° for 2 at larger periods, even 
between OUT and SUK. There is a definite phase shift 
between OTJ and EIS, but because the short period 
variations are suppressed southeast of EIS a clear 
meversallwsSonot Seems Im Stack 3 there is asphaseushitt 
between NOK and LON. This fact and a too large separa- 
tion between magnetometer stations may be the major 


reasons why the conductor is less obvious in Stacks 2-4. 


4.3 Maps of Fourier spectral components 


Maps of Fourier transform amplitudes and phases 
were introduced in the first array study by the Alberta 
and Dallas groups (Gough, 1973) and serve two main 
purposes. Firstly they give a great deal of qualitative 
and some first order quantitative information, and 
secondly, they form the basis of any further quantitative 
interpretation. 

Four sequences were selected for Fourier 
transformation from the three data sets. The event of 
January 16 was transformed over the time interval 14.00- 
18400 8GeMe ho and the -event of wanuWary 2-22 over the 
interval 11.30-19.3N G.M.T. on January 21. The event of 
January 28, which was the one most successfully recorded, 
was transformed over the intervals 14.30-22.30 G.M.T. and 


17.39-22.30 G.M.T. Figures 4.6 and 4.7 show amplitude 
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Fourier amplitude spectra in the period 
range 20-250 min for 14.00-16200° GINO 
January 16, 1972, for the stations 
OKASwEl Ss ands Ee 


Figure 4.6 
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Figure 4.7 Fourier amplitude snectra in the period range 
20<240 min for 17.30-22230° GeMUIe y danuaryecoe 
1972... for the stations OKA. DEI wandananh: 
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Spect%a fOratnes time periods 1400-18-00. Gator. 
January 16 and 14.30-22.30 G.M.T. January 28 for three 
Sitdcions. 

Periods with spectral peaks for Z and an 
associa tedspeak inseither X or Y, butivery often in both, 
were chosen and contour maps of amplitudes and phases 
at these periods were prepared according to techniques 
described by Reitzel et al. (1970) and Camfield et al. 
(1971) eandiealso discussed in Section 2.2. Eleven cets 
of maps, each containing three amplitude and three phase 
maps, were drawn. The period range covered is 21-171 min. 
The polarization of the horizontal field was elliptical 
for the 11 periods mapped. The azimuth of the major 
axis ranged from N 60°W to N 17°E. 

The sets of Fourier maps shown in Figs 4.8-4.13 
have been chosen as representative of the 11 sets of 
maps. These maps cover the period range 21-120 min and 
an azimuth range of N 60°W - N 13°E. Only maps transformed 
from the two sections of the January 28 event are presented 
because the stations in the eastern part of the array 
recorded this event most efficiently. In the western part 
of the array where the three events were equally well 
recorded, the contour maps show the same features in all 
11 sets of maps, but in the east the maps for the January 


28 event show more detail because of the availability of 


more data. 
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Figure 4.8 
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Fourier transform amplitudes (nT/c/min x 10) 
and phases (min) at period 21 min for 14.30- 
Zand Gollalba. Wanliary> 2ere 972 | VES timate: eto 
the standard deviations in the X, Y and Z 
components are respectively 30, 30, 16 nT/c/min 
fom the ampli tudes and 55, 3.5 eandi2. bainiteror 
tie ohases. lhe volarpzation of vine nortzontal 
fielders indicated: 


JANUARY 28, 1972 T=28 MIN. 


X AMPLITUDE x X PHASE 


Y AMPLITUDE Y PHASE ae 


Z PHASE 


Figure 4.9 Fourier transform amplitudes (nT/c/min x 10) 
and phases (min) at period 28 min for 17.30- 
Oca OirtGeM edt. J AnUary con 972.) ES Caliates: at OT: 
the standard deviations in the X, Y and Z 
components are respectively 25, 29 and 13 nT/c/min 
fOmethe samo latudesseand) Wro.. loand 275 "nin “tor 
the phases. The horizontal field polarization 
is indicated. 
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JANUARY 28, 1972.T=35 MIN. 


bourier transform amo latitudes: (nd /c/min 2 10) 
and phases (min) at period 35 min for 17.30- 
C20 Gan ales, anrtiary 2c JA9722 2 Estimates © for 
the standard deviations in the X, Y and Z 
components are respectively 24, 24, 16 nT/c/min 
for ‘the amplitudes: and, 5 mandi nine tor 

the phases. The horizontal field polarization 
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Figure 4.11 Fourter transform amplitudes (nT/c/min x 100) 
and phases (min) at period 54 min for 17.30- 
Pew Ganal., vanltanry co.) 9/2. SEst inmates 7oy 
the standard deviations in the X, Y and Z 
components are respectively 26, 30 and 14 
hige/mueror che anpiitudes=and 1." | owand elas 
min for the phases. The horizontal field 
polarization is indicated. 
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JANUARY 28, 1972. T= 79 min. 
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Fiaqure 4.12 Fourier transform amplitudes (nT/c/min x 100) 
and phases (min) at period 79 min for 14.30- 
Bos G4 Me | ewan Ua tveccow 197 26 we ESita mate Saanor 
theses tandard devigtivons in the x. Yeard 2 
components are respectively 32, 37 and 20 
nT/c/min for the amplitudes and 1, 2 and 2 min 
for the phases. The polarization of the 
horizontal field is indicated. 
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JANUARY 28, 1972. T=120 MIN. 


Fiaqure 4.13 
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Fourier transform amplitudes (nT/c/min x 100) 
and phases (min) at, period 120° min for 4. 30- 
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All the X amplitude contour maps show a 
prominent anomaly in X elongated in an east-west direction. 
Me stations SUK; EIS and SEH are "situated near the 
maximum of this component. 

Seven of the Y component contour maps have a 
maximum at SUK and EIS coinciding with the maximum in X. 
The horizontal field was polarized north-west - south- 
Gast (ergs 4090. 4510, 45112 4913) stor most of “these 
periods. This component increases to a local maximum at 
the south-eastern edge of the array (Figs 4.8, 4.9, 4.11 
and 4.12) in five of the ¥ amolitude contour maps. For the 
two largest periods mapped, T = 120 and 171 min, the Y 
component amplitude maps show a minimum at KWA in the 
eastern part of the array (Fig. 4.13). Both these periods 
have the horizontal field polarized slightly west of north 
(Fig. 4-13). Maps for five other periods with a polariza- 
tion east of north show a local maximum at KWA (Fig. 4.12) 
in the Y component. 

The Z amplitude contour maps are more complicated 
than those for the X and Y components. A maximum in Z 
Occurs at the stations SUK and OUT on the western edge of 
the array for all the Z amplitude contour maps compiled 
(Figs 4.8-13). A second smaller maximum in Z is present 


at REH near the south-western corner of the array. 
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The Z amplitudes in the eastern half of the 
array are invariably smaller than those in the western 
half. In the eastern half the Z amplitude anomalies are 
also polarization dependent. With the horizontal field 
polarized west of north SEP, NOK and HAP show a local 
maximum and SEH and KWA minima as in Figs 4.9-11 and 
4.13. These anomalies become less prominent for the 
horizontal field polarized more towards the north-east 
as in Figs 4.8 and 4.12. The maximum at SEP and NOK 
for the horizontal field polarized north-west - south-east 
is in agreement with a curved conductor trending east-west 
in the western part of the array and curving north-east - 
south-westward in the eastern part of the array (Gough 
et al., 1972). Some features in the eastern part of the 
array may not be very well defined due to the instrument 
Spacing. 

The phase maps for the X and Y components appear 
to represent mainly the normal field and to be uninfor- 
mative of internal structure. The Z phase maps are, 
however, of considerable interest. Although the Z phase 
values could have large errors in the eastern third of 
the array because the Z amplitudes are small there, all 
the Z phase maps and Table 4.2 show that there are 
large phase changes between the stations OUT and SUK, 


TSU and GOB, NOK and HAP, and between SEP and GWE. In 
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Table 4.2 


Phase Differences for Z 


T min OUT-SUK TSU-GOB NOK-HAP SEP-GWE 


Ad deg Ad deg Ad deg Ad deg 

January 16 Ze) an ete +144 uel FAO) 
Bl -162 - 94 - 69 

January 21 73 - 43 - 77 - 88 
114 -156 -123 -113 

171 _-160 te + 8 

January 28 21 +169 oy +215 nail Res: 
28 +154 +139 oe Te 

35 +144 AG +169 she SiS 

54 -173 -150 -154 =113 

79 -159 -142 -172 - 46 


120 = 159. chs! -139 -139 
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terms ote thessimple assumption of a current flowing 
under the maximum gradient in Z phase and along the 
contours of the maps, Figs 4.8-13 and Table 4.2 suggest 
a current along a curved conductor passing between OUT 
amad=SUK near EIS; SEH and KWA; and thus substantiate 

the suggestions offered by the magnetograms and amplitude 
Mapsee Gough, Lilley and McElhinny (1972) nointed out 
that a curved conductor produces Z amplitude anomalies 
which are strongly dependent on the azimuth of the 
polarization. Since actual substorm fields are 
elliptically polarized, the phase difference between 

the X and Y components will enter into the problem 
(Bennett and Lilley, 1972; Gough, McElhinny and Lilley, 
1974). Due to these complications it is necessary to 
work in terms of a parameter of Earth response to 
incident fields which preserves phase and amplitude 
information. Transfer functions between the components 
meet the need. Transfer functions as used by Schmucker 
(1964, 1970a) were thus computed in order to clarify the 
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The array was situated at a magnetic latitude 
at which the normal Z component in substorm fields is 


relatively small. If it is assumed that the Z variation 
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fieldseaneemainlyeanomalous., 1t 16. 4Usti tiable to cal= 
culate transfer functions from estimated normal horizontal 
components xy and Vi to Z. Although this method avoids 
inversion of the full transfer function matrix which 
relates each anomalous component to each normal component 
(a procedure which requires a full separation of normal 
and anomalous fields), it should be handled with care 
because a systematic Z-H correlation in the normal fields 
could affect the calculations adversely (Schmucker, 1964). 
The transfer functions between the Z component 
(assumed anomalous) at a particular station and the 
estimated normal X component were calculated using the 
RewatlTons 701 eG. cocy.. AS d1S.cCUuSssedsin section 2.0. the 
real part of the transfer functions Zy and Zy were 
combined to yield an in-phase ‘induction vector’. This 
induction vector was as usual reversed to be in the same 
sense as a Parkinson arrow (Parkinson, 1959, 1962) so 
thaueartepoints toWards a good conductor in which the 
current flows in phase with the normal field. Similarly 
the imaginary parts of 2, and Zy were combined to yield 
a quadrature-phase ‘induction vector’. To be consistent 
Withethe in-phase induction Vector, the sense of the 
quadrature-phase vectors were reversed as well to point 
towards conductors in which the current leads the normal 
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phase induction vectors as shown in Fig. 3.11 were not 
reversed. 

Nine 120 min events which had been recorded 
in three components at 13 stations were selected from 
one three main events studied. These are Wisted in 
Fables 43. Various combinations of six of the mine 
events had been recorded at 18 stations. Transfer 
functions computed from nine and six events recorded 
Ate cers stations yielded Simidarenesuits a5 lhe 
induction vectors derived from all nine events, except 
at RUN, EIS, NOK and GOB where induction vectors are 
computed from the events of January 16 and 28 and at OTd 
where they are computed from the events of January 16 and 
2leare oresented in Fig. 4.14. Four different periods are 
shown. For each spectral term of each event an estimate 
of the normal field in X and Y was made by taking the 
mean cosine and sine Fourier coefficients over all these 
Sit dit i-Onis 

The reversal in both the in-phase and quadrature- 
phase induction vectors clearly indicates a linear curved 
conductor underlying the array. The quadrature-phase 
induction vectors are generally smaller than the in-phase 
induction vectors. Along the western margin the quadrature- 
Hidseninduceion vectors for | = 43omin and T= Ize min 
tend to indicate the east-west conductive structure, while 


the in-phase induction vectors show a strong deflection 
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Table 4.3 


Two-hour Events Used in Calculating Transfer Functions 


Date GoM t. 


dl Uta View Onna We Outee 14.00 = 16700 


MSU MUA 72 gs Tie eS Oe ale 3.0 


January 2o, 1972 oes Oe an 130 
Meee SO. =o os 0. 


"ee et 


yey? ie? yreemes 


ot rans 


Figure 4.14 In-phase (solid) and quadrature-phase (broken) 
induction vectors representing the transfer 
functions Zy and Zy at periods 14 "min, 26 mine 
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towards the coast. This deflection of the in-phase 
induction vectors is to a lesser degree also present for 
io= 14 min and. |= 26.mins but the directions of the 
quadrature-phase induction vectors are more scattered 
for these periods. 

This separation of the in-phase and quadrature- 
phase induction effects seems to indicate a structure 
with high reactance off the western edge of the array, 
while the east-west structure in the western half of 
the array has a low inductive reactance compared to 


its resistance. 


4.5 Normalized anomalous fields 


Profiles of normalized anomalous fields were 
constructed for the section A-B along the western margin 
of the array (Fig. 4.1) and the section C-D near the 
eastern edge of the array. Since anomalous fields cover 
most of the array, the stations OKA, NAM, MAR, ARA were 
arbitrarily designated as representing normal field 
components. The mean X amplitude from these stations 
was used to define the assumed inducing field along the 
profile AB and, as residuals, the anomalous horizontal 
components x. at other stations along the section. For 
the profile CD the horizontal component 45° west of 


north of the same four stations was designated as 
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representinge the inormal’ or inducing field. his 
component is parallel to the induction vectors in the 
area. The mean Z amplitude at these stations was 
deducted from Z at every other station on the profiles 
to give an estimate of the anomalous vertical field Za: 
Under these assumptions Lai Ay and Na were determined 
along the profile. Since the amplitudes of the anomalies 
Showed a marked dependence on the polarization of the 
horizontal inducing field, Fig. 4.15 shows the 

profiles along AB for which the inducing field was 
polarized west of north and Fig. 4.16 shows those for 
which the inducing field was polarized east of north. 

The most prominent feature in these profiles 
is the reversal in Z./X, between OUT and SUK and the 
associated maximum Of; ee (unfortunately the X 
component at OUT did not function). 

It is clear that the anomaly has larger ampli- 
tudes especially in Zi! for the horizontal field 
polarized with an azimuth west of north, which is in 
agreement with a northward swing in the conductor 
further east. The ratio of Nal is not very different 
in the two diagrams, except for T = 35 min where ee 
is very large. This is most probably because the azimuth 


Gfeits horizontal tield polarization 1s the furthest west. 
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Figure 4.15 Normalized anomalous variation field components 
| i Pie and X aie along section AB (Fig. 4.1) 


for Be ieee 1 min (event of January 16, 1972) 
28, 35%and 120 (event of January 26. 1972). 
For all these periods the horizontal field 
Doli¢rvuza tone 1s west ofanoy tar 
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The amplitudes of es and Z,/%, decrease 
with increasing period, with T = 35 min again an 
exception. 

For the section CD the average horizontal 
field component in a direction 45° west of north for 
the stations OKA, NAM, MAR and ARA (the "normal" stations 
for profile AB) was designated as representing the 
"normal" or inducing field, He This component of 
the horizontal field parallels the induction vectors 
along the section. The average vertical component at 
the above-mentioned stations represents the "normal" 
vertical component Zs These values of ue and ZL, were 
subtracted from H (the horizontal component 45° west 
of north) and Z at every station along the profile to 
give estimates of ne and Zs: Figure 4.17 shows eight 
normalized anomalous profiles nis and raya at eight 
different periods together with the average polarization 
Of the Norizontal freld nor that period. sihewproniies 
are arranged so that those that have polarizations west 
of north are on the left side of the diagram and those 
with polarizations east of north on the right. 

All profiles show a maximum in Zo/H, ES) Ole 
near the station NOK, with a decrease in Zaft further 
South=-eastward.. Except Tor 1 = 54 min, all profiles 


indicate a maximum in Us to the south-east of the 
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Normalized anomalous variation field components 
oe and Heyy along the sections CD (rig, 4.1) 


at periods 23 min, 51 min (event of January 16, 
1972) = 24 mins 2oomMins. soamin. o4:°ninee so oni n 
and 120 min (event of January 28, 1972). The 
polarization of the horizontal field 1s shown 
for each period. 
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station MAK. The main difference between the profiles 
Of. di frerentshorizontal’ field polarizations 1s" the 
maximum in reals atenOK for all profilesewith polariza- 
tions west of north. At a period of 79 min with the 
horizontal field polarized NNE there is an indication 
of this feature, but for all other periods having 
inducing field polarized east of north, there is no 
maximum in HE Ale at NOK. This behaviour most probably 
indicates that the north-east trending branch of the 


conductor consists of more than one parallel conductor. 
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The magnetograms, maps of Fourier spectral 
components, transfer functions and normalized anomalous 
fields all show a conductor with axial line approximately 
cheap OS itd Of) IMGT Cated lI acd 4c Ne DOS ION son 
the conductor is not very well defined near the eastern 
edge of the array. When the attenuation in Z at LON, 
the anomaly pattern along profile CD and the high Y 
amplitude values near the south-eastern margin of the 
array for certain periods and polarizations are 
considered it 1S guite possible that there are more than 
one conductor in the eastern half of the array. This 
Can only be established by further studies im tne area. 


The east-west trending conductor joins the NE-SW branch 
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which most probably continues as far south as LON, near 
SEH. Although the polarization sensitive behaviour of 
the anomaly precludes meaningful two-dimensional modelling 
a maximum depth can be estimated from the half-width of 
the anomaly defined as the width of the ete anomaly at 
half-maximum value. A line current at depth d produces a 
X field of width 2d at the half-maximum value, anda Z 
field with extrema 2d apart. The actual current cannot 
be deeper than the line current best fitting the anomaly 
widths just defined (Gough, 1973p). | 

In this case the width at half-maximum value 
Of ae on profile AB is 250 km and the Z extrema are 
100 km apart. This difference is due to the fact that 
the X component at OUT was not recorded. This reduced 
the definition of the an anomaly and the station 
Spacing of the array was not designed to detect such 
narrow anomalous zones. These half-widths will indicate 
maximum depths of between 125 and 50 km. I am inclined 
topaccept a maximum depth closer to the smaller depth 
especially if the separation of the in-phase and 
q@uadrature-phase induction vectors is considered. This 
separation of the induction effects seems to indicate 
thavethemeast=west conductor ts ataaslitthospierrvcedentn. 
The conductive structure joins the African rift system at 


the eastern point of the array where it underlies the 
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Middle Zambezi Rift Valley containing Lake Kariba. This 
rift zone is the south-westward extension of the 
Luangwa Rift. Reeves (1972) found that the Kalahari 
seismic axis is indeed a well aligned extension of 
Zambia's Karroo-filled Luangwa Rift and the induction 
anomaly coincides with the seismic active areas. In 
Botswana the linear features in the Okavango lie parallel 
tomtnes conductor, DUt in South-West Africa the conduc- 
tive structure cuts obliquely across the general grain 
in the Damara Geosyncline. De Beer, Gough and Van Zijl 
(1975) reviewed the available geophysical and geological 
data in the area and concluded that the conductive zone 
marks an extension of the African Rift system along old 
weak zones in the lithosphere. 

The deflection of the in-phase induction vectors 
towards the west could be due to the proximity of the 
continental edge and the Atlantic Ocean, but could also 
be caused by anomalous structures associated with the 
prominent north-south faults of Karroo to Recent age 
indicated in Fig. 4.2. Korn and Martin (1951) showed that 
most of these faults are Still active, If the in=phase 
response 1s considered, however, 1 1s clear that it 1s va 
structure with a high inductive reactance and is most 
probably a subcrustal conductor associated with the 


continental edge. 


ee E> aw 

stat. .palign! San: taba eat ley t34% Feedee a 
salt te. wo! OO SPRS: byeee come Puce ato zi vhs 

“Nia whe fed? GrOPRARTeL? tevesk ee oe pire) 

2G Sika eo Sosaet of eTKs stusted 7 

=i p37 eared Bart tt-oo wes vhatanek | 


M4 
7 


» (meat Wee aot tontes ¢(eRONE 
Aik =} fas 2 Sebemaee seh, oe pase? of. 
tens ei Yee .1etcuenes on?- 63 : 


45 22025 esugneriee oid 


' moet Te ae si:anypeee exe of at 
4 , deo aie) Vewerven SRtGT) 
Meodate -he?) be bi sanen King ano Ge AP SPED. 
: , apbonsia’s fie 2am 7 
Fr es vn? 2enps i56h, . 
teach iasenae it) ‘ 29 noryoer tay ae | ; 
sf Sivoo Seem Sf) 2Sepwas _ 
wens of Pusey sia Ge eens totmagtina | 7 
i> Oy tart lise ou RAK “sl betes: vt, 
ij of Suvset ju eat arabe drag ident ae 
Ae een reO( ‘syt4c@ Pat roan Gage ple 7 bes a 7 
seyde—nt B37 sy ‘tad)! Gere yi cs 


rr. d4 Seed AAS 2f ot il 


Pig 27 one sonbs 256" weal a 


7 a - 


aie ag Fw hadiR don tau p Bier iT ; 


161 


The work in this chapter has been the subject 
of a paper by De Beer et al. (in press). 
The results of this array study will be 


discussed in more detail wm Chapter 6, 
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CHAPTER V 


PS OAUDYNOR Sint ELECTRIC CONDUC TIVINY eS TRU CIURE 
IN THE KAAPVAAL CRATON 


5.1 Electromagnetic induction in a plane-layered Earth 


The dimensions of magnetometer arrays in mid- 
latitudes are usually smaller than the spatial wave- 
lengths of the inducing magnetic fields. When this is 
the case and there is no indication of major lateral 
inhomogeneities in the area under investigation, the Earth 
can be approximated by a plane-layered model and the 
problem becomes that of solving the normal equation (1.14) 
for plane layers. The solution of this problem is quite 
~ general and has been treated by Price (1950, 1962), by 
Wait (1953) and by Schmucker (1970a), among others. 

In this induction problem the Earth is represented 
by a model with N-1 conductive layers above a conductive 
halfspace (Fig. 5.1). The model occupies the region 
z > 0 in a right-handed Cartesian coordinate system. The 
J2thetaver Clot elseos Noi N) OG threkhess k = 2541-25 
has conductivity oe: The air above the surface is layer 
9 and oe is assumed zero. A periodic primary field with 
exp (iwt) time dependence is assumed. 

Price (1962) showed that the solution of the 
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E = exp(iwt)Z(z) F(x,y) (Seat) 
where F(x,y) = (so. 2 ae 0) (Soca 


2 2 
~ 4 = ke = 10 (5.3) 
X vy 


The parameter k is a measure of the non-uniformity 
of the primary field, and is given by k = 2n/X where A is a 
linear dimension of the source field. Thus k = 0 implies 
a uniform source. The function of depth Z(z) satisfies 


the equation 


2 
Sas =o 4 jwuo)z (sua) 
dz 
2 
d*z _ 12 
or welt Kae (ey) 
7h 
with Ke = Ke + iw 9 (Geaea 


Demiseecleanetiiataeds +s) hase sOlUTTONG BroOrst Ope. tnes norm 
Sink). ocos( kx). exp (ik), and sin(k x + ky) where 

These solutions involve single values of 

k which may be regarded as individual modes of the source 


field. 
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The magnetic field is obtained from the 


equation 


eee Shoah 


and is given by 


foe] 


Z ap 


2 F 
- 37 by? k=ZPy exp liwt ) Coy 7) 
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In the j-th layer the solution of eq. 5.5 is 


Z5(k,t,z) = A, (k,t)exp(-K,z) + B i(k, t)exp(K,z) 


(523) 
dz ,(k,t,z) 
d ——>+—_————— = - K.{A.(k -K. - B. : 
an ra 7 4 ,t)exp ( K 52) B exp (K 5z)} 
(5.9) 
g. dz, 
Sie. Gi ast Zs = - K ame 
A exp(-K.z) + B exp(K;z) 
Bere g 52) a AseXxp Kz = B exp K 42 10) 


, 2 ] 

Eneetias case g; (2) = G tz) of Schmucker (1970a, 
ela. WSR a 9; (2) = ae (alice of Weaver (1973, p. 275) and 
apy 4) Pe g, as used by -SChnMUCKe ta Gyles pwaO4 

At any interface the tangential magnetic and 


electric field intensities are continuous. This implies 
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an Zs(Z.44) ; 4254142544) 
Zz dz 
Dut Talso that 
94) 4ia1) _ Sger' 2547? (Seales 
K jt] 


From eq. 5.10 we have at the lower boundary of the 
j-th layer, i.e. interface jtl, 
A. =Kee7 omy bie eZee 
GP AiZe gece ince 52541) je eee 
may ict 2 a ° a . . e 
Ji aia A exp ( Kee B exp(K,z. 17) 


aq abe 


(513) 


The coefficients A. and BS can thus be eliminated using 
eqs 5.10 and 5.13. The positive and negative arguments 


— 


are combined to hyperbolic functions to yield 


95(2544) + tanh K (2544-2) 


ea c 
g2\2) = ~ (5a 1:4") 
1 [Ss q;(Z5,)) tanh K (2544 Zs) 
with 2. aace 244] 


Equation 5.14 corresponds to eq. 82 of Weaver 


iT] 


(HiS7 So) Settina: z Zs in this equation gives 95 (25) 


J 
in terms of 95 (25,1). Equation 5.12 is used to obtain 
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the basic recurrence formula for pnlane-layered conductors 


K 595 4q'2 541) + K 541 tanh K, hs 


Cie 2) iate Kari aslo metaniin aah ie (ort oy) 
Apr Keay + Ri94 47 2 54] tan iF j 
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with h; = 254] ~ 2 = the thickness of the j-th layer. 
The above equation is the same as the expression for oF 
Onepem 1 O4 of Schmucker (1971). 

The incident field disappears at great depth 


(Zeteyiso that the term B. in the Noth layer (the 


N 
halfspace) must be zero and 


Jy (Zp) Ss (CSral6y) 


The value of q, (9) can thus be obtained by 
starting at the surface of the halfspace where Gy (Zy) = ] 
and calculating successively the q;(z,) by means of 
Gidie 5b. )o% 

Two limiting cases can be considered. Firstly 


when the skin-depth in the j-th layer is large in 


comparison to the wavelength of the source field, then 
Kaen (yeners) 


and the incident field penetrates through tnis layer as 


if it were nonconducting. 


The second case arises when the wavelength of 
the incident field is large in comparison with the skin 


depth of, say, the j-th layer. Pyonmeed. 655,60 
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K; = VIOUS; = ¥ (5.18) 


with 6 the skin depth defined for the j-th layer as 


Shee ae (5.19) 


the attenuation of the field in this layer is governed 
only by the skin depth and is independent of k. 
Im terms of eas, 2520, eoe2- wand) 527 0:9 


& ap 
a = exn(iwt) {A ,exp(-K;z) + B exp(K5z) bay 


; - oP 
and Tout = exp (iwt){A,exp(-K.z) = B exp (K 32) } 57 


Gray 


so that the magnetotelluric impedance (eq. 1.24) at the 


earth's surface is given by 


E iwu.g, (0) 
ie on Oa all (5.22) 
H K 
y ] 


and is independent of P. 


From eq. 5.7 we have that at the earth's surface 


H 
oe tet E (5.23) 
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but this relation is dependent on:P. 


Furthermore we have from eq. 5.3 that 


One om (a 
Roya lee <a Goma) 
OX oy 


oH aie aaeG (5.25) 
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again at the earth's surface. This equation gives a 
relation between the vertical variation field and the 
horizontal field gradient and is independent of P. 

Schmucker (1970b, 1971) introduced the 
inductive scale length C (w,k) - the response function 
used by Weidelt (1972) - and defined it as 


C(w,k) - " (5.26) 


Since we are working at the surface the 
SsubSermipt of G is dropped in the further discussion, Let 
will imply C,. 

The magnetotelluric impedance (eq. 5.22) reduces 


to 


he iwu,C(wsk) (Sea 72) 


(AS. 2) 
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andseq. 5.25 can be written as 


H 
oy Sqr = C(w,k) (o528) 
potent SS) SINS cn 
OX oy 


Schmucker (1970a) and Kuckes (1973a) independently 
introduced a simple inversion technique for one-dimensional 
problems. This method yields a first estimate of the 
conductivity at depth from the response function at the 
Surface. These authors assume a uniform half-space of 
conductivity 5 underlies a layer with a low conductivity 
and thickness h. In this model the incident field will 
penetrate the top layer with negligible attenuation for 
frequencies of interest in array studies. The penetration 
depth in the half-space will be small compared to the 
spatial wavelength of the incident field. 

From the limiting cases discussed in connection 


with eqs 5.17 and 5.18 we have for this case that 
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where 55 is the skin depth in the second layer. 
Equation 5.15 yields for a two-layer ‘case 
Katze eteken adil) Dunant 
Perle cee Zi ] (5.30) 
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and when the limiting value for Ko 1s tmeroduced 
g, (0) = kC (55329) 
6 6 
2 re ys 
where G = h + Paina ues (5533) 
a complex length. Since k = K, 
a, (9) 
C = "4 = Ca, k) (sno) 
] 
the resvnonse function of eq. 5.26. 
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Thesgeneral strategy of the inversion) 1s to determine sc 
at the surface by means of equations such as 5.27 and 
Seco Inesequations 5.35 are tnen Used to get for each 
frequency component a depth and conductivity of a uniform 
"substitute conductor". This depth and conductivity are 


meaninaful only when the argument of C(w,k) for a specific 
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frequency lies between-7 and Q, yielding a positive 
depth h. For such a frequency a multi-layered structure 
is indistinguishable from and therefore replaceable by a 
uniform conductor at depth h. An underlying assumption 
is that (h + 5.) is small compared to the wavelength of 


the source field, 


* op 
With hh ~~ Cosy 
the depth of a perfect substitute conductor 
* 
ip SRG) S72) 


and a modified avparent resistivity can be defined as 


eerie ie 2wu {Im(C) } (5.38) 


5.2 Considerations for interpreting observed data 


The theory developed in Section 5.1 applies to 
Specie sources fleld configurations -and speciiticc 
geological conditions. In practice the investigator 
is faced by the problem that he seldom knows the source 
field configuration and actually sittudies the conductivity 
structure to learn more about the geology. Under certain 
conditions, assumptions can be made regarding these 
parameters to facilitate the computation of theoretical 


models. 
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Lilley (in press) considers the effect of 
different elementary P functions on eq. 5.23. For a 


travelling wave one can take P of the form 


P(x) = exp(ikx) (5.39) 
Wirtne eh real. Then 


P Pay 
dP/dx sa ca (5.40) 


and eq. 5.23 reduces to 


= —— = ikC(w,k) (S544) 


This equation has been the basis of what Lilley terms 
traditional geomagnetic depth sounding, as developed for 
example by Price (1962), Whitham (1963), Caner et al. 
(1967), Cochrane and Hyndman (1970) and Caner (1971). 
NAften the modulus of both sides of the equation is taken, 
SHUS MUOSHNGetie 1 dla 1tsestanit cance sm inmtnl sequal on 
a specific Knowledge of k iS required, With P of tne 
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or for the general case with P a function of both x and y 
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Wath the total horizontal field. This dives. an 
estimate Of kK. Tne modula of both sides of this equation 
ane or ten Cdken, again losing the Significance of tne 


i-factor. For a travelling wave it is however valid to 
dH 


replace the ik “factor in eq. 5.41 ‘by alt, to get 


2 ee (5.44) 


For standing waves one can take P of the form 


P = sin kx, when 


Pp Se otan kx 


odes nk (5.45) 
so that eq. 5.23 becomes 

H 

7 = -kC(w,k) tan kx (5.46) 

x 


This means that for standing waves eq. 5.23 is no longer 
independent of x. The use of eq. 5.23 to interpret 
observed ly ratios is thus valid only for a single 
mode travelling-wave source field. 

The magnetotelluric and array-gradient equations 
(5522 and 5.25) seem not to depend on the question whether 
the source-field is of travelling-wave or standing-wave 
nature. In these cases one has to bear in mind that tne 


relevant field components will vanish at nodal points. 
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The value of k does not enter in these equations in the 
Same €xplicit way aS in eq. 5.23, but. does enter in the 
estimation of the K, values. This point will be 
considered further on, in the discussion. 

Schmucker (1970b, 1973) discusses the general 
induction problem as seen from a particular surface point 


Le randae es 


ii 3 
The first scale length, Li; refers to the spatial con- 


in terms of three relative scale lengths L 


figuration of the source field and can be thought of as 
the spatial non-uniformity of the overhead source field. 
This can be for instance the height of an overhead line- 
current jet or the reciprocal of the wave number k for a 
Sinusoidal field. The second scale length, Lo» is related 
to the degree of lateral uniformity of conducting matter 
beneath the surface. Schmucker (1973) assumes that the 
halfspace z > 0 is subdivided at vertical boundaries 

into different layered conductive structures. The length 


L, is defined as the distance between the point of 


2 
observation and the nearest vertical boundary. Weidelt 
(1972) showed that im the special case Of Sutficientiy 
uniform source fields the positive real part of the 
inductive scale C(w,k) can be interpreted as the weighted 
mean depth of the induced currents in a stratified half- 


space. The limiting value for k + 0 is a convenient 


measure for the depth of penetration at a given frequency. 
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Schmucker (1973) thus defined its modulus for the 
relevant stratified section as the third scale length 
such that L; = |[C(w,0)|. Obviously L, increases with 
increasing period. 

If we assume a quasi-uniform source field, 
that the observation points are far removed from lateral 
boundaries and that the subsurface approximates a 
1? Ly So L3. This 
is a one-dimensional induction problem because, as 


layered model, we have the case where L 


Schmucker (1973) points out, for the area under considera- 
tion neither the exact configuration of the source field 
nor the conductivity in adjacent stratified sections 
matters. For this case one can set k = 0 to obtain the 
field relations for the so-called Caaqniard case from eqs 
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Equation 5.49 can be used to obtain an estimate of Co G:) 


from the observed vertical field and horizontal field gra- 
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dients at d certain frequency. Tt is worth noting that 
H 


PrOMEOGR HO .23 7 = 0 when k = 0. To circumvent this 
me 

problem Schmucker (1973) makes allowance for a small 

Sinusoidal modulation of the source field with the 


constraint that 


ke 


C(w,0)| = p= << 1 (5.50) 


Using this he then defines a locally effective "wave 
number" (personal communication between Schmucker and 


Eilleyes hove) 
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This means that both sides of eq. 5.49 are normalized 
WhLiMreSpecy Ou | Fl eand SOudth fers Tome une suse sOt 6d o-4. I, 
where the observed eae ratio is compared to the computed 
Valles sikclw.k).. Schmucker considers the reciprocal 


of Ke doo auscele length for the overhead Source. 
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A study of the properties of C(w,0) as derived 
promreq.) 5-49 scan be helotul in examining the valadiwey of 
the assumption of a nlane layered earth with induction 


by a quasi-uniform source field. Firstly C(w,0) has to 


dRe[C(w,0) ] 
dW 


the depth of induced currents changes smoothly with 


vary smoothly with frequency and < 0 because 
frequency, with the lower frequencies having larger 
penetration depths for a given conductivity distribution. 
EROMPeds =5..53 1h Ss cleakethat RepC(a.0 ) 12-90 and 
Imfe(w.0) =< 05 1f Cw 0) as determined from the obsierved 
data does not have these basic properties, some or all 

the assumptions regarding a one-dimensional model are 
violated. Another consequence of a departure from these 
properties will be that the Schmucker-Kuckes inversion 


technique cannot be applied. 
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compared to a global model 
The short period limits for studies of the 

radial distribution of electrical conductivity in the 
Earth and the resulting poor determination of the 
conductivity structure for depths less than 400 km have 
been discussed in Section I.2.1. “Ine northern part oF 
the 1971 array was situated on the Kaapvaal craton and 
the magnetograms, maps of Fourier spectral components 


and induction arrows discussed in Chapter 3 indicate that, 
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of all the areas studied, this area seems to be the least 
affected by lateral inhomogeneities. It would thus seem 
that these array data offer an opportunity to determine a 
conductivity-depth profile for an Archean shield area. The 
results of the deep Schlumberger soundings in the area 
(Section 1.4) indicate that the geological structure does 
indeed approximate a plane-layered situation. It should, 
however, also be borne in mind that even at the north- 
east of the array H(t) bears a strong resemblance to 
H(t), suaqgesting that much of the ie variation field is 
an anomalous field resulting from induction by HY in the 
east-west conductor to the south of the array. 

In principle the best approach to the problem 
of obtaining a conductivity model for the Kaapvaal craton 
should be a direct inversion technique such as those 
described by Bailey (1970) and Weidelt (1972). These 
two methods, however, require precise data over a broad 
frequency range, while in this case only data with a 
limited precision over a very limited frequency band are 
available. Weidelt himself advises: “it appears that the 
best way to handle geomagnetic induction data is still to 
interpret them by a set of homogeneous layers and to 
introduce, if necessary, further preconceived model 
assumotions" (Weidelt, 1972, p. 285). The modelling in 


this study was therefore performed with the aid of the 
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Conductivity models investigated in one- 
dimensional model studies. Model 1 is the 
"standard" model constructed from conductivity 
models by Banks (1972) and Van Zijl and Joubert 
(1975). Model 2 is derived from Model 1 by 
increasing layer thicknesses between 250 and 600 
km by 50%. Model 3 is derived from Model 2 by 
decreasing conductivities between 225 and 900 km 
DY a factor. on. OF 
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theory developed by Price (1962) and Schmucker (1970a) 
andedtscussed=invSection 521, 

Because of the effects of the Cape Fold Belt 
conductive structure on the data and the assumptions 
involved in the modelling techniques, it is clear that 
the present data allow, at best, a comparison with 
existing models and precludes any meaningful independent 
modelling. A good model for comparison is therefore 
required. Despite the limitations in the data of Banks 
(969s 1972) discussed in Sectionel<2.1% his=models for 
the radial conductivity structure of the Earth are 
considered to be the best available. The conductivity 
models derived from deep Schlumberger soundings are the 
best available for the crust and uppermost mantle 
(depths < 200 km) in the craton (Figure 1.3), where 
lateral variations of conductivity can be reduced to 
acceptableslimitss (his modell@and thatvas Banks (1972) 
were therefore combined to form the "standard" model for 
comparison and is depicted as Model 1 in Figure 5.2. 

The most severe problem in this type of modelling 
TSPto "ge tean estimate of the spatial wavelength of the 
source field. Price (1962) estimated that for geophysical 
applications involving sources of natural origin the 
spatial wave number would lie in the range 1.5/7 x 107! m7! 


POW Aes 1072 m7 |. Banks (1969) states that for mid- 
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latitude stations away from the auroral and equatorial 
electrojets generally k < 2/Re =o | ax 102m |. For k in 
this range the approximation C(w,k) ~ C(w,9) holds 
(Larsen, 1973; Schmucker, 1973). This means that the 
first scale length Ly Pos EWithe theses cation CHURMihtiG ee owl 
as observation point, the nearest major lateral inhomogeneity 
observed in the array area (under the Cape Folded Belt) 
is about 600 km away. There may be, and most probably 
are, closer vertical boundaries outside the array. On 
the present data, however, Lo TS aboutio0@ skims 

The third scale length Ly is defined as 
Mod[C(w,0)]. It will be shown that this quantity varies 
from about 409 km at a period of 20 minutes to about 
1,000 km at a’ period of 300 minutes. This means that 
at best we have the situation that Ly >> Lo, La and do not 
actually have a one-dimensional case. It is however 
instructive to treat the problem as if Ly; Lo >> L. and 
to observe the effects of the lateral inhomogeneities. 

As a first study the classical geomagnetic 
depth sounding approach was investigated. The observed 
a ratio was modelled using eq. 5.41 and the algorithm 
described in Section 5.1. The value of k in eq. 5.41 was 
determined by means of eq. 5.43 using only the moduli of 
H and H.. This approach thus assumes a plane-wave source 


x 
configuration. Straight lines were fitted by least squares 
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to the HY Values dt the stations PEI, KOOP. DEA eSCH. ICHRs 


dH 

VRY (Fig. 3.1) to determine —~ and to the H_ values at 
aX sH y 

OLI, KUR, CHR to determine ae . The estimates of 


the standard errors of the slopes were combined with the 
Cnrorseditier “CoO, get. the standard errorssit ke vA wel canted 
average of thesk estimates at 10) periods: (3248) 76)) x 1077 m7 
was used in the calculations. This assumes that the k 
values for different periods belong to the same population 
and may not be a valid assumption, but certainly is better 
than guessing values. It will be noted that this value of 
k coincides with the Banks (1969) upper limit of 
3.1 x 107/m7!. A11 stations on the Kaapvaal shield as 
outlined in Fig. 3.1 were used to obtain an average 
a rato. VA rom Fig. 523, 1 te iseclearethat ons, quanti ty 
does not vary smoothly with period, an indication that 
the assumptions do not held) forall the data used. ine 
figure also shows the |=" values for the standard model 
(Model 1 in Fig. 5.2). At short periods the model 
approximates the observed ratios, but has too small a 
slope to match the observed ratios at the longer periods. 
To determine the effects of perturbations to 
the model the thickness of the layers between 250 and 600 
km depth in the standard model was changed to 150 km to 


Optaimemodel 2 (fig. 5.2) and the cellculation yepeated, 
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Figure 5.3 Observed le values for stations on the 


Kaapvaal shield in relation to computed 
values for. Model 1 and Model 2. The values 
of k are for the one standard deviation 
Limits: 
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Ine results are shown in Fig. 5.3. The slones of the 
model curves are <¢til] too small. Model 3 in Fig. 5.2 
was constructed by increasing the resistivities of the 
layers in the depth range 225 to 900 km a factor of 10 
relative to Model 2. This model fits the data better 
than the standard model (Fig. 5.4). 

The effect of the crustal model used was 
investigated by extending the lower crustal - upper mantle 
resistivity (1,300 ohm.m) right to the surface. The 


effect for k = 3.4 x 107! ! 


is shown by the broken line 
CURVED Tn (Figs, 5.4 “and Tt ts “obVI0USs that near=surface 
conductive layers have only a slight effect at the 
Shorter periods and none at the longer periods. 

The normalized anomalous profiles in Figs ~ 3. )2 
and 3.13 indicate that the fields at tne sotithern’ stations 
on the shield (to the north-east of PRI) contain sizable 
anomalous parts due to the conductor to the south. The 
Srrectl on including only the noreiernmos stations in 
the computation of an average 2 was assessed by calcu- 
lacinaethisuratio: for CHReand VRY. sai Nel res UItGi ss simon 
Timp dees 5) ine relation to theoretically aluves Tor stne 
standard model and Model 2. Now the standard model fits 
the data for T < 70 min quite well and Model 2 fits all 
theedata satisfactorily. [tt 1s interesting to note 


that in the range T < 70 min the approximation that 
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Figure 5.4 


H 
Observed |=“ values for stations on the Kaapvaal 


shield in relation to computed values for Model 3. 
The broken line curve indicates the effect of 
extending the lower crustal-upper mantle 
Vesa S Cid LY. one). U0 Sonne m to the surface of the 
Earth for k = 3.4 x 10-7m™ 
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Figure 5.5 Observed Fa ratios at the stations CHR and 


VRY in relation to computed values for Model | 
and Model 2. 
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Ly» Lo >> L3 is more appropriate than for longer periods, 
because L3 increases with increasing period. The above 
calculations were repeated using the approximation of 
eq. 5.52, but with k as determined from the moduli of 
tievhield components, The results are almost identical 
to those discussed above. The biggest uncertainty in 
thisetype Of analysis is the estimation of k,” If for 
instance the largest k (6.2 x 107/m7!) obtained in the 
computation of the average k is used in the computation, 
the standard model gives a good fit to the data. 

A method that avoids the estimation of k is: 
therefore needed. The array gradient method is such a 
method. The observed response function was computed for 
different periods using eq. 5.49. In this case ae H 


oH 
and Hy were treated as complex quantities and aa was 
determined along a S-N line through KOF, DEA, SCH, CHR, 
OH 
and eVRY {( Fig. 3.1) and sy along an W-E line through 


Ria hURGHR.| (ihe waverage He value at the stations 


X 


VEN, VRY and CHR in the north-eastern corner of the array 
wasused in the computation of C(w,0). The values as 
computed at single periods were found to contain large 
standard deviations, partly because of the small Hes 
Values, but mostly because of the errors*in the gradients. 
A frequency-band averaging technique was therefore devised 


to decrease the errors in C(w,0). In this process response 
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functions were evaluated over a band containing up to 15 
Fourier spectral estimates at periods near 20 min and two 
estimates at 300 min and then averaged over the band. 
This procedure is appropriate because the response 
function chances very slowly with frequency and the 
C(w,0) within a limited frequency band can be considered 
to be an estimate of the same quantity. A similar 
technique is used in magnetotelluric studies in the 
estimation of the tensor impedance elements (Word et al., 
970s Hermance; i973). 

The bands were selected to include portions of 
the spectra with high energy in Ho LOuINcreasea the 
Signal to noise ratio. The bands were also chosen over 
a period range where the polarization angle varied no 
more than about 30 degrees. Stability in the polarization 
direction was assumed to indicate stability in the 
source field. Changes in the polarization direction 
within a band may be valuable to eliminate the effects of 
local inhomogeneities, but for the present study stability 
in the source field was considered more imoortant. A 
reason for this was that the C(w,0) values varied a large 
amount within a band if the associated horizontal field 
polarization direction changed by a large amount. The 
Fourier coefficients represent averages over a time 
interval and changes in polarization with period may 


indicate contributions from more than one source. 
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The frequency-band averaged values of Mod[C(w,0)] 
are shown in Fig. 5.6. The response functions are plotted 
at the geometric mean period of each band. Only the 
moduli are plotted because in only one band were Re[C(w,0)] 
positive and Im[C(w,0)] negative as they should be (Weidelt, 
1972). This agrees with the conclusion, from the earlier 
analysis in terms of scale lengths, that the structure 1s 
not in fact one-dimensional. The penetration depth increases 
with increaSing period as it should. The theoretical 
response UNG tGIOnSs. Tor Model Mimand Modelme3=of F 1Geeo.c 
are superposed on the data and the latter clearly fits 
the data better. This is consistent with the results 
from the study of 2. 

In some of the bands investigated but not shown 
in Fiq. 5.5 the values of C(w,0) were unstable because 
slopes with about equal magnitudes, but opposite signs 
were added in the denominator of eq. 5.49. [his problem 
can be avoided by rotating the time domain data to 
Components parallel and perpendicular to the long axis 
of the polarization ellipse. This will mean that the 
data have to be Fourier analysed to obtain the polariza- 
tion parameters before the rotation can be done. The 
rotated time dcmain data are then Fourier analysed and 
the gradients cbtained in the new reference frame. In 


the case of a linear or near-linear polarization only 
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Figure 5.6 Frequency-band averaged values of Mod[C(w,0) ] 
in relation to computed values for Model 1] 
and Model 3. The vertical bars indicate the 
standard errors and the horizontal bars 
indicate the width of the frequency bands. 
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one derivative has to be taken which will enhance the 
accuracy of C(w,0). In this study there are too few 
stations in the study area to carry out this process. 

Although the data used in this analysis have 
some properties that indicate a deviation from a one- 
dimensional structure, the results suggest that the 
conductivity profile in the approximate depth range 
200-800 km under the Kaapvaal shield differs from that 
of the Banks model. The observed fields require a slower 
rise im conductivity with depth. Magnetometer array 
Studies are well suited for studies in this depth range 
which is deeper than that for Schlumberger and 
magnetotelluric studies and is near the short period 
limit of radial conductivity studies. A magnetometer 
array located on the Kaapvaal and Rhodesian shields might 
Helowto establish an electric conductiviry profile ior 
Chases nield region. 

With more stations removed from major induction 
anomalies, a trend surface could be fitted to the 
horizontal field components which would lead to better 
estimates of the gradients and also C(w,0). A serious 
attempt should also be made to develop a Backus and 


Gilbert type inversion technique for this type of problem. 
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CHAPTER VI 


THE CONDUCTIVITY ANOMALIES IN SOUTHERN AFRICA IN RELATION 
TO LOCAL GEOLOGY AND GEOPHYSICAL PARAMETERS AND THEIR 
POSSIBERSOCAUSES 


6.1 The Southern Cape induction anomaly 


The orinel pal result of the 1971 array “study 
discussed in Chapter 3 was the discovery of the east-west 
trending conductor under the west central Cave Fold Belt 
andsaeepepartao; the Karroo basin (Fig. so. 14). Inds 
conductive structure was so far from being exnected that 
most of it passes to the south of the southern corner of 
the triangular array. Nevertheless it is clear that a 
large highly-conductive structure elongated east-west 
underlies the deep Karroo basin and the Cape Fold Belt. 
The phase relations and frequency response of transfer 
functions discussed in Chapter 3 indicate that the 
conductor is most likely in the lower lithosphere or 
asthenosphere. The relations of this conductivity anomaly 
to other aeophvysical parameters in the area and to the 
local geology are very interesting and warrant a short 
discussion. 


6.1.1 Heat flow data 


Conductivity anomalies in the crust and uppermost 


mantle can be related to temperature, as in the western 
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Figure 6.1] 
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Published heat flow data (Gough, 1963; Carte and 
Van Rooyen, 1969; Chapman and Pollack, 1975; Von 
Herzen and Vacquier, 1967) for southern Africa in 
re Llatinan tio, tectonic nvavinces of 4 a 1 2.8l Thea 
labelled tectonic provinces are in chronological 
order: A the Kaapvaal shield, B the Limpopo Belt, 
C the Irumide Fold Belt, D the Namaqualand-Natal 
Belt, E Damaran-Katangan age province and F the 
Cape Fold Belt. The star indicates an average of 
20 heat flow measurements in Lake Malawi (Von 
Herzen and Vacquier, 1967) and the other eight 
values north of 20°S are the recently reported 
values for Zambia (Chapman and Pollack, 1975). 
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U.S.A. (Gough, 1974), or to composition (Garland, 1975). 
Figure 6.1 shows some published heat flow data (Gough, 1963; 
Carte and Van Rooyen, 1969; Chapman and Pollack, 1975; Von 
Herzen and Vacquier, 1967) for southern Africa Superimposed 
on the tectonic provinces of Fig. 1.2. There are no heat 
flow data from the Cape Fold Belt, but at least three 

heat flow measurements on the conductivity anomaly just 

to the north of the fold belt and one further east just 

to the north of the folding (Gough, 1963). These four 

heat flow values have a mean of 60 mWm~° 

GieHEUe== | aecal Chee = Ato mW) which is 1 ae more 
than that for the other six values on the Namaqualand-Natal 


Belt. There seems to be no higher heat flow over the 


anomaly. 
Gelee  eGraVvaty data 


A large negative isostatic gravity anomaly exists 
Inesouth-eastern South Africa dt theseas tern end .o7 sthe 
Gane Fold Belt (Fig. 6.2)... lt reaches =a minamum of -800%g-4- 
(meq yes tls Re Once ns, mgal) of which about) 2000g.u; 
can be ascribed to Cretaceous sediments (De Beer et al., 
1974). This anomaly is elongated along an east-west Bocuse 
In repvorting this anomaly Hales and Gough (1960, 1961) 


exnlained it as being caused by a relict root or thickening 


of the crust which had provided Airy type compensation for 
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Magnetic, gravity and conductivity anomalies in 
southernmost Africa in relation to the Cape Fold 
Belt. The maximum and minima of the Beattie 
magnetic anomalies shown in Fig. 6.3 are 
indicated. The aravity anomalies are departures 
from an Airy model of isostatic compensation 

with crustal thickness 30 km. The inset map 

is adapted from Hales and Gough (1960) and the 
gravity units are in mga] (1 mgal = 10 g.u.). 

The arrows represent the negative of the in-phase 
induction vectors for a period of 128 min for tm 
five stations nearest to the anomaly. The 
approximate axial line of the conductive 
structure is indicated. 
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a mountain range since eroded away. It was pointed out 
that a mountain range which would just be supported by 
the root would be 30 km wide and 1.4 km high. Since 
large stresses are necessarily associated with isostatic 
anomalies of this large magnitude and area, a sequence of 
failures was predicted by Hales and Gough (1960) leading 
Finally vosup lift near the “anomaly axis and reduction: of 
the isostatic anomaly and the stresses. There are no 
large scale isostatic anomalies in the western half of 
the Cape Fold Belt which could be explained by postulating 
that the process of unlift is more advanced there than in 
the east. Sthnge (1935) showed that large vertical 


movement had in fact occurred in the western region. 


Oe iece neromagne tic data 


The most prominent linear geophysical feature 
in South Africa is the Beattie Ridge magnetic anomaly 
(Beattie, 1909) which lies to the north of and, for most 
of 1ts 900 km length. parallel to «the Cape Fold Belt 
(Kidseesceand 6,35). 8 hen Malm anomaly coms) sus) OT ta 
maximum with ampolitude ranging from about 100 to 500 nT 
(foteloantensity) tlanked over large distances by long 
linear minima. Inspection of the 1:1,000,000 Bouguer 
anomaly map of South Africa (Smit, Hales and Gough, 1962) 


reveals no associated gravity anomaly. Deep Schlumberger 
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soundings across the anomalous magnetic structure (Van Zijl, 
Bahnemann and De Beer, unpublished observations) suggest 
that the causative body is in the basement beneath the 
Karroo sediments. Model studies by F. Bahnemann (personal 
communication) indicate that the anomalous magnetic structure 
is most likely in the middle or lower crust. Du Plessis 

and Simpson (1974) found that at their eastern extremeties 
the linear east-west trending Beattie anomalies are 
truncated by the continental slope, These authors showed 
from simple depth calculations that on the continental 

Shelf the bodies causing the anomalies lie at depths 

varying between 0.6 km and 3.0 km. In this depth range 
seismic refraction studies have been interpreted as 
indicating Cape-Karroo rocks (Ludwig et al., 1968). This 

is not consistent with the results from the deep Schlumberger 
soundings further west. Du Plessis and Simpson (1974) 
believe that the causative bodies are late-Karroo dolerite 
intrusions rather than Precambrian intrusives and metamor- 
Hes. WS Sioleiielstell mo tee Como wake Milaee Wolk 
situation further west in the array area, because there 

the Beattie Ridaqe anomaly lies well south of the southern 
Karroo dolerite occurrence. There are no Karroo dolerites 
intruded into the western part of the fold belt and the 
whole fold belt actually shows a total absence of any 


sianificant igneous or regional metamorphic activity. From 


See I27 Oo SA. On. aac mony, Ea 


Pies See ees Dee 

respi: i" Fewaaedy > a 

ont ‘dthhhhasangress? Seeme, ap ebod vShnind on Jad 

reaqag, ovenanted 7 We evihpie [etch ,2tnemtbse IOVS: 
ce “TISAPGM | eo.) Exp ne ofa feds ul §o Fen (1oh¥0> Faun uN 

svt ct ©. keows cweedT the Sibbie 2\t ni xtodtl Feo . ” 
teixs / shee" Tla97) VEER aaa (nper? noeamte bbe 

34 af lentne )atovecn me! uae te tigw-Jene thantt id i 


One 217 ft] 4 : 23RngaG ney oe aa sit ye gti: ! : j 
ink a-Si NS +2 16 ies kes tat aTqate mort 
+ gee py ae Be (usa (eotbad lied tote oa 


dress oth WL “ede eee ll a0 psewied dues 
524 sinh cbioete Meliseyiey ee. =) 

F . pases ore seed antes ton 
a jnstel choo tonet 

() wear pe baa. By eat Youll. ~srey vsidTet rondo aloe 
| awe ghey By eatin ets. ‘aca aneet toa 

bes Beer ink bat meee wea ay | fc + 


‘ ay ‘ 
ha Pome b= res J412 > SA 28 ae bpaiag ts ea bY 
vi = 


tits aad a die meee pat gt hae mht 


4 


. 


(eAlai Say BOA HP Ae ol | hs cane mabe salah af tha ey a 
‘ oe ie | i 
i TOV VEO Sa Fee Oe 
mane 
S72 tie yeae oi yee sre, vee Siz es 


a es wingees ea ie 
vs eer ae titan irom 


206 


Figs 6.2 and 6.3 it is clear that the Beattie anomalies 


do not continue under the fold belt. 


Oe lea the Cape Fold Belt and its origin 


ihe Cape Fold Belt consists of two folded 
chains, a smaller north-south belt along the western 
edae of the fold mountains and the principal east-west 
chain along the Cape Fold Belt marked in Fig. 6.2. The 
folded sediments ranqe in age from Late Ordovician (Cocks 
et al., 1970; Cramer et al., 1974) to late Permian 
(Haughton, 1969) but the main phase of folding was in 
the Triassic (Haughton, 1969) and involved the east-west 
chain. The less severe north-south trending Cedarberg 
folding occurred in pre-Dwyka times (i.e. before the 
Middle-Carboniferous). The lithology and palaeontology of 
the Table Mountain and Bokkeyeld Groups (the two lower 
stratigraphic units in the Cape-Karroo sedimentary sequence) 
argue for deposition on a continental slope with marine 
conditions to the south (Rust, 1973). The first signs of 
a southern land mass (Ryan, 1967) and pre-Dwyka folding in 
the east-west fold belt (Haughton, 1969) are observable in 
Wittebera (Lower Carboniferous) strata. There is now good 
evidence that this southern land mass was not simply an 
uplifted are of the continental craton as was widely 


accepted (e.g. Haughton, 1969; Newton, 197.3). o.bute an a sd.anid 
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arc type of orogenic welt possibly related to a subduction 
zone with associated trench south of the present coast- 
line. Near the southern edge of the Karroo basin gray - 
wackes rich in volcanic fragments and tuff layers are 

found in sediments that were transported from the south 
(Martini, 1974; Elliot and Watts, 1974). Since there was 
no volcanic activity in the Cape Fold Belt, the volcanoes 
must have been situated seaward of the present south coast, 
in Gondwanaland of Permian times. 

The origin of the folding in the Cape Fold Belt 
is still a matter of controversy. Some geologists (e.g. 
Newton, 1973, 1974a,b) argue that the fold belt does not 
resemble a typical "Plate Tectonics" orogen in respect of 
structural style, sedimentation or igneous and metamorphic 
activity. Newton (1973) therefore proposed a gravity- 
folding model controlled by therstructure of the: pre- 

Cape rocks. 

Other geologists (e.g. De Swardt et al., 1974) 
suggest that the Cape orogeny is of Alpine type generally 
ascribed to gravity sliding away from a median zone of 
excessive uplift, with movement of the cover over a more 
rigid basement, which is accompanied by folding and 
thrusting. It is further pointed out by these authors 
that the best explanation for orogenies of Alpine type 


is that of collision of continental lithospheric plates. 
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An imoortant factor in the assessment of the 
two proposals is the relative position of southernmost 
Africa in the pre-drift Gondwanaland. Geophysical evidence 
concerning the south-eastern continental margin of South 
Africa (Francheteau and Le Pichon. 19725) Seruc tons and 
DugMLessa se 197 32 scruttoni. 197ca,bee0UuLPLlessi cmand 
Simpson, 1974) as well as drilling results from the Falkland 
Plateau (preliminary reports from Scientists aboard Glomar 
Challenger for Leg 36 of DSDP, 1974) strongly suaqgest an 
Upper Triassic-Lower Jurassic Gondwana model in which the 
Falkland Plateau as part of the South American plate 
bordered the African plate and extended as far as Durban 
along the present coast (Scrutton, 1973a). The Gondwana 
reconstruction would then be like that of Smith and Hallam 
(1970) with the Antarctic Peninsula perhaps positioned 
further south. Patagonia and the Antarctic Peninsula 
would be sicuated» to thersouth of Africa, taking directions 
with reference to present-day Africa. In his book "Our 
Wandering Continents" Du Toit (1937) described this fold 
belt crossing southern Gondwanaland, and it is clear that 
the Cape Fold Belt represents only the northern edge of 
this Gondwanide Orogen. 

The available evidence seems to support the 
hypothesis that during late Palaeozoic-early Mesozoic 


times on oceanic plate underthrust Gondwanaland along its 
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then southern boundary. Finally a continent-continent or 
more likely continent-island arc collision occurred 
producing mountain building which lasted from the Permian 
to middle Triassic. In South Africa the most obvious 
result of this collision is the Cape fold mountains. 
Over-folding to the north is a prominent feature in these 
mountains indicating that the pressure was directed mainly 
from the south. 

ihe v,absence of volcanic activity in the Cape 
Fold Belt is consistent with a subduction zone with trench 
off the present southern coastline of Africa. The geology 
and nalaeontology of the Falkland Islands and results of 
deep sea drilling on the Falkland Plateau (preliminary 
report from the Scientists aboard the Glomar Challenger for 
Leg 36 of DSDP, 1974) indicate that the islands were also to 
the north of such a subduction zone. 

Padula et al. (1967) presented geological evidence 
for a subduction zone (they refer to an eugeosyncline) that 
bordered the western rim of present-day Patagonia, possibly 
from Devonian to Triassic times. The likelihood of this 
subduction is well supported by the presence of a paired 
metamorphic belt extending along the Pact Vee coasts or SOUTH 
America from Santiago de Chile to Tierra del Fuego. The 
high-pressure metamorphic zone is situated on the western 


side of this belt and the age of the metamorphism is 
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Permian (Gonzalez-Bonorino et al., 1970). In Patagonia 
late Palaeozoic to earliest Mesozoic magmatic ACTIV TTY. 
both intrusive and extrusive, has been demonstrated for 
the North Patagonian Massif and southward to the Deseado 
(essijeialperne 1068s halpernveteg lea 1972) oe Rubudvune 
strontium whole-rock isochrons and mineral dates indicate 
igneous activity from Carboniferous to earliest Triassic 
for the South American sector Gondwanide Orogen. 

In West Antarctica there exists similar evidence 
for subduction. The deformation is of early Mesozoic 
age in the Ellsworth Mountains and through West Antarctica 
there are late Palaeozoic and early Mesozoic intrusions, 
which suggest that the Gondwanide Orogen crossed this 
part of the continent (Halpern, 1968; Dalziel and Elliot, 
1973). The upward increase in air-fall tuff beds in 
Triassic strata of the Beardmore region suggests that 
the centres of volcanism were becoming more vigorous or 
slowly encroaching on the Beardmore area (Elliot and 
Wattise 1974). A subdictton to the southe of Permo-iriassic 
Gondwanaland is thus firmly established. 

A foldina model based on continental (or 
continent-island-arc) collision offers an explanation for 
the neaative isostatic anomaly associated with the Cape 


Fold Belt, because the roots of collision mountain belts 
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This model, however, does not offer a direct 
explanation to the presence of the Beattie magnetic 
anomalies or the induction anomaly which is actually 


of primary importance in this discussion. 


Oto) Proposed causes of induction anomaly 


Gough (1973c) suqgested that it is possible 
that the high conductivity zone is associated with 
anomalously high temperatures and partial melting and 
that the conductive body could in fact be a linear 
region of ascending mantle material - a ridge-like 
Structure on the more conductive zone in the mantle. 
The isostatic anomaly may then be the aravitational 
Signature of this same east-west trending region of 
high temperature in the upper mantle. If this is the 
case further work should indicate that the isostatic 
and induction anomalies coincide in the eastern part of 
the Cane Fold Belt. The heat flow data in Fig. 6.1 
neither support nor deny this suggestion. In the 
earlier discussion it was noted that the main phase of 
folding occurred during Triassic times. The folding 
itself therefore will be unrelated to any possible 
thermal anomaly in the mantle because it is hardly 
possible that such an anomaly Would persist for tatime 
of order 10° yr. Recent uplift in the region, however, 


could be related to a present thermal structure. 
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Despite this Newton (1974a) suggested that 
Gough's proposal supported his model for the formation of 
the fold belt. Newton proposed that such linear mantle 
bulges would offer a driving force for the vertical 
movements of the fault-defined blocks which his model 
requires. He also suggested that the conductivity 
anomaly should lie to the north of the isostatic 
anomaly because the isostatic anomaly lies on the 
depressed southern side of the major fault bounding the 
Cretaceous-filled Algoa basin. This is in contradiction 
to Gough's explanation of the isostatic anomaly and 
would also imply that the thermal anomaly persisted for 
about oe vi. 

A set of data that has not been used in this 
discussion is the radiometric ages for the basement rocks. 
Nicolaysen and Burger (1965) pointed out that a south- 
westerly traverse starting on the Kaapvaal craton would 
pass from granitic gneisses older than 2,600 Myr old 
gneisses of the Namaqualand-Natal Belt and then through 
the Malmesbury aeosyncline intruded by 600 Myr old 
granites and capped by the Palaeozoic-Mesozoic Cape Fold 
Belt. In this extensive age study of the Namaqualand- 
Natal Belt they suggested that the thermal event that 
caused the metamorphism could have been a linear incursion 


of heated mantle below the Namaqualand-Natal metamorphic 
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zone. They further suggested that the continent grew by 
the accretion of younger mobile belts like the Malmesbury 
geosyncline. Put into modern plate tectonic language 
this could mean that a plate was underthrust beneath the 
Gondwanaland of about 1,000 Myr ago along its southern 
boundary to produce the high temperature belt that caused 
the 1,090 Myr old metamorphism. This process also 
accreted material onto the continent to extend the 
continent southward and in this process formed the 
Malmesbury geosynclinal sequence into which the 600 Myr 
old Cape Granites are intruded. A similar, but much more 
extensive process has been suggested for the formation of 
most of Western Canada west of the Rocky Mountains (Monger 
Sito Ween Ou 2 Jr 

The Beattie magnetic anomalies are in agreement 
with such a linear accreting process. From Figs 6.1 and 
6.3 it is clear that these anomalies roughly parallel the 
southern boundary of the Kaapvaal craton. The aeromagnetic 
maps of British Columbia also show linear anomalies of 
comparable amplitude parallel to the coast (Haines et al., 
1971; Stacey, 1973). The conductivity anomaly therefore 
could be caused by the transition from the highly metamorphic 
crust and associated uppermost mantle to the much younger 
(at least, 1,600 Myr younger) crust and uppermost mantle 


of the Malmesbury geosyncline. In its exposures to the south 
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of the fold mountains the Malmesbury rocks consists of 
avis: shales, graywackes, limestones, basic lavas and 
tuffs, quartzites, phyllites, schists and with higher 

qrade metamorphic rocks in the areas intruded by granites. 
In general the metamorphism was low-grade (Haughton, 1969). 
Highly magnetic peridotite incorporated into the crust 
could give rise to the Beattie magnetic anomalies. This 
model would explain the association of the induction 
anomaly and the static magnetic anomaly. From Fig. 3.14 

it is clear that the Karroo sediments are thin to the 

north of the station BEA (Beaufort-West). which as on the 
northern edge of the induction anomaly, while to the south 
of it the Karroo basin deepens dramatically and the folding 
begins. This is most probably due to the difference in 
mechanical properties between the older lithosphere 

north and the younger lithosphere south of Beaufort-West. 
In this case the northern edage of the conductor would 

mark the southern boundary of the Kalahari craton (Clifford, 


1967). 


Coulee Conc husa On 


This discussion is necessarily speculative 
because some of the crucial data to determine which of 
these hypotheses, if any, is the correct one are still 


lacking. The required data will include an extension of 
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the array coverage further east, heat-flow studies in 

the Cape Fold Belt and seismic crustal studies to determine 
if such a drastic transition from the Namaqualand-Natal 
Belt to the Malmesbury geosyncline do exist. Geological 
investigation of the problem is severely hampered by 

the thick cover of Palaeozoic-Mesozoic sediments in the 
important areas. The vroblems concerning the Cape Fold 
Belt and its relation to the basement geology and 
geophysical parameters in the area form part of the South 
African contribution to the Geodynamics Project. 

On the present data, which are in many respects 
incomplete, the writer prefers the hypothesis that the 
induction anomaly is due to a compositional difference 
between the lithosphere in the Kalahari craton and the 
younger lithosphere further south. This hypothesis 
seems to me more consistent with the available data than 
the proposal of Gough (1973c), but a decision is impossible 
at this stage. 

In Chanter 3 a NW-SE trending anomaly was 
described, seen at the stations OLI, PRI and PET for 
horizontal fields with north-westerly polarizations. This 
anomaly coincides with the Doornberg-Brakbos fault system 
which marks the boundary between the Kaapvaal craton and 
the Namaqualand Belt (Viljoen, 1975) and which ends in the 
North Cape Fold Belt. This fold Belt shows up in the gravity 
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maps for South Africa (Smit, Hales and Gough, 1962) and 
Botswana (Reeves and Hutchins, 1975) and runs north 

from the stations UPI and OLI in Fig. 3.1 to just south 
of the Okavango Delta in Botswana (Fig. 6.4). In Chapter 
4 evidence was given for a conductor running east-west 
through South-West Africa and turning northeast in 
Botswana to enter the mid-Zambezi valley. It was shown 
that induced currents were laraqer for NW-SE polarization 
of the horizontal field than for NE-SW polarization. The 
North Cape Fold Belt may link this conductive zone to the 
Doornberg-Brakbos fault zone, in a structure generally 
consistent with the observed polarization dependence of 


the induction anomaly. 


6.2 The conductive zone in South-West Africa and 


Botswana 


The most Significant result of tne 19/2 earray, 
study in South-West Africa, Botswana and western Rhodesia 
is the discoverv of a zone of concentrated induced electric 
currents indicated in Figs 6.4 and 6.5. The station 
spacing defines the boundaries of the conductive zone 
within about half a degree, except at its western 
extremity where they are determined within a quarter degree. 


The work in this section has been reported on 


by De Beer, Gough and Van Zijl (19 753) 
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Fiqure 6.4 


The conductor in relation to seismicity and 
tectonics of southern central Africa. Small 
dots ,.°2)s5o Mex. 4e Ss lharidewd OC seein omen SH cal 
Squares intensities < 33 large squares > 3. 
Faults (bold lines) are from Geological Map 
of Botswana (1973), Gough and Gough (1970), 
Geological Map of South West Africa (1963), 
Provisional Geological Map of Rhodesia (1971) 
and International Tectonic Map of Africa (1968). 
Geological boundaries are from the geological 
maps of Botswana and South-West Africa and 
Haughton (1969). The bathymetric depths are 
in metres. 
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Fidqure 6.5 African ritts ‘and sive lidcmarneicel ationen omcne 
conductivity anomaly (compiled from a map by 
McConnell (1972) and other sources). 
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The half-widths of normalised anomalous fields 
and the induction vectors indicate that the conductor is 

at crustal depth. From east to west the conductive zone 
continues the south-westward line of the Luangwa-middle 
Zambezi Rift (Fig. 6.5) and bends south of the Okavango 
Delta to run nearly westward across South-West Africa. 

In the north-eastern part of South-West Africa and all 

of Botswana the Kalahari sands of Tertiary or younger 

age cover the solid yock and hide almost all tectonic 
features. The Kalahari basin was formed in Tertiary times 
when large parts of southern Africa went through phases of 
intermittent and differential uplift intercalated with 
periods of peneplanation or pediplanation (Haughton, 1969). 
Boreholes reveal Mesozoic Karroo sediments or volcanics 
under much of the Kalahari basin indicated in Fig. 6.4 
(Reeves, 1972; Reeves and Hutchins, 1975) but have not so 
far penetrated the Kalahari sediments in the Okavango 

Delta (D. Hutchins, personal communication). Major linear 
features with a north-east trend are found near the 
Okavanao and have been interpreted as faults (Reeves, 

1972: Geological Map of Botswana, 1973; Scholz et al., 

in press). These lineaments continue the line of the 
middle Zambezi-Luangwa Rift, an important linear structural 
element in the tectonic framework of central southern Africa. 


The middle Zambezi valley containing Lake Kariba, is an 
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asymmetric rift preserving Karroo rocks downfaulted along 
about Middle Jurassic faults (G. Bond, personal communica- 
tion) between Precambrian shield rocks. Although the 
fault pattern is reasonably well known (Cox, 1970; Gough 
and Gough, 1970) little is known about the mineralization 
and water bearing qualities of the faults (G. Bond, 
personal communication). In the Luangwa Rift we also 
have Upper Karroo sediments and volcanics downfaulted 
against basement rocks, with the age of faulting generally 
the same (Drysdall et al., 1972). Several igneous centres 
are known along the Luangwa Rift (Cox, 1970), including 
Cretaceous carbonatite volcanoes (Vail, 1968; Bailey, 
1961). There are also carbonatites and kimberlites in 
the middle Zambezi (Lee, 1974) and a cluster of diamond- 
bearing kimberlites near Orapa in Botswana (Fig. 6.4). 
As early as 1921 Gregory considered the Luangwa-middle 
Zambezi Rift as an extension of the East African rift 
system better developed and known in Kenya and Ethiopia. 
Others, like Du Toit (1954) and Bailey (1961) used 
geological arguments to propose south-westward continuation 
of the Luanqwa-middle Zambezi Rift into Botswana. 

To test the hypothesis that the conductive zone 
indicates a continuation of the known rift system of 
Africa other geophysical parameters in the better studied 


rifts further north are compared to the same parameters in 


the study area. 
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6.2.1 Seismological evidence 


The two northern branches of the East African 
Rift system (Fig. 6.5) are characterized by continuous 
belts of normal faulting and graben structures and the 
occurrence of shallow earthquakes (e.g. Baker et al, 1972; 
Wohlenberg, 1975). The Eastern and Western Rifts, however, 
differ from each other in several other respects (Maasha 
and Motnar, 1972). “The seismic activity of the Eastern 
Rift is much less than that of the Western Rift (Gouin, 
1970; Fairhead and Girdler, 1971) and increases to the 
south. Molnar et al. (1970) suggested that much of the 
deformation in the Eastern Rift occurs aseismically. In 
contrast to mid-oceanic ridges and rift structures seismic 
activity extends to some 100 km beyond the main rift 
flanks (Wohlenberg, 1975; Bath. 1975) 

Considerable seismic activity with magnitudes 
up to 6.7 occurred in the Okavango Delta region in the 
period 1949-1955 with 90% of the events in late 1952-early 
Wo5s (Cane and Oliver, 1953: Oliver, 1956). “Reeves (1972) 
mapped 50 earthquakes in Botswana over the period 1965-1971 
and found a concentration of epicentres in the Okavango 
Delta and scattered events in the central Kalahari. Reeves 
believes the broad Ghanzi Ridge of post-Damara sediments 
separates the two seismic regions. In the Lake Kariba 


region a low initial level of seismicity rose several orders 
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of magnitude as the lake filled, culminating in three 
events of magnitudes near 6 in 1963 just after the lake 
level peaked. Triggering of renewed activity on the 
existing faults, by incremental solid stress rather than 
water pressure, was suggested by Gough and Gough (1970). 
SMice 190s "activity Was econtinued with very slow, if any, 
decline, in a manner consistent with the hypothesis that 
the region is now seismic in a steady state. Mechanisms 
for two of the large events show normal faulting with one 
nodal plane striking N 25°E parallel to regional faults 
(Sykes, 1967; Fairhead and Girdler, 1971). These and 
-other fault plane solutions for earthquakes in eastern 
and southern Africa (Banghar and Sykes, 1969; Fairhead 
and Girdler, 1971: Maasha and Molnar, 1972) suggest 
tensile deviatoric stress in a WNW-ESE direction. Fairhead 
and Girdler (1971) showed that in the period 1963-1970 
earthquakes in southern central Africa were located in a 
broad seismic zone running southward from Lake Tanganyika 
tomamsooint just west of Lake Kariba. Except im tie Lake 
Kariba area the middle Zambezi-Luangwa Rift was virtually 
aseismic through 1963-1970. 

Epicentres from several sources (Fernandez, 1972; 
Fernandez and Guzman, 1973; Rhodesia Meteorological Services, 
1973; L. M. Fernandez and J. A. Guzman, unpublished) for 


the period 1965-1973 have been pO tted iit kgs eGc4n emunenny 
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distribution supports the hypothesis that the Luangqwa- 
middle Zambezi Rift continues to the south-west into 
Botswana. The focal depths of the earthquakes range down 
to 50 km (Fernandez, 1972; Fernandez and Guzman, 1973). 

A micro-earthquake study was done in the Okavango 
region from mid-September to mid-December, 1974 (Scholz 
etvalas in) press). © The activity was largely confined to 
a narrow band running northeast along the southeastern 
edge of the Okavango Delta and further north-east towards 
the Zambian border. The rate of seismicity was found to be 
about the same as that in the main rift valleys of East 
Africa. A composite focal mechanism for well located 
events near the southeastern boundary of the delta 
indicates normal faulting on north-east striking planes. 
One nodal plane dipping 60° to the northwest strikes 
Daralilel tol the fawlts- in) the area” Scholz et als, (in 
pressye suggest that! the delta is at the tip of a Zone: of 
IncimpTenteritting that tol lows the older” Nid=Zambezi 
Luangwa Rift. 

West of the Okavango the conductive zone runs 
a few degrees north of due west to the exposed part of the 
Damara Geosyncline in South-West Africa (Fig. 6.4). The 
zone of high conductivity cuts obliquely across the 
WSW-ENE grain (Haughton, 1969; Clifford, 1967) of the 


high-arade metamorphic rocks of the eugeosynclinal Swakop 
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Facies of Upper Proterozoic-Lower Palaeozoic ages 
(Clifford, 1967). The miogeosynclinal OUTIO Factas= to 
the north has generally E-W strikes in its eastern 
exposures and gradually assumes a more northerly strike 
westward. The conductor runs roughly parallel to the 
folding in the eastern outcrops of Outjo Facies rocks 
and the boundary between the eugeosyncline and the 
MUOaTOsy ne line. tt also parallels the northern limit 
of the Damara regional metamorphism (Clifford, 1967). 
No recent major faults are known to run Nena lel fe ime 
At least four earthquake epnicentres were located close 
to the northern edge of tne conductor for the period 
1965-1973 (hiv 6.4). Ani association of the ‘conductor 
with the boundary between two stable crustal blocks 
seems possible. 

Seismological evidence for the mantle under 
thesbast African rifts shows low P velocities and attenua- 
tion of Sy north of a discontinuity near latitude 4°S 
(Fairhead and Girdler, 1971; Maasha and Molnar, 1972; 
Baker et al., 1972; Wohlenberg, 1975). For a path 
Addis Ababa-Nairobi, Rayleigh wave phase velocities are 
low for a continental path (Knopoff and Schlue, 1972) and 
resemble those in the Basin and Range province of western 
North America. Elsewhere in Africa phase and group 


velocities of surface waves and travel times of anes S, 
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and Ly indicate that the structure beneath most of Africa 
is typical of stable continental reqvons. (bloch, Guede. 
1969; Gumper and Pomeroy, 1970). 

Bloch, Hales and Landisman (1969) studied 
Rayleigh-wave phase velocities in southern Africa and 
found lower phase velocities west of a line through 
Pretoria and Bulawayo than along that line and east of it. 
In the western region which includes the magnetometer 
array area the Rayleigh-wave phase velocities are very 
Similar to those in the central United States and eastern 
Prrce on na which are both non-shield continental regions. 
This is consistent with the tectonic pattern of southern 


Agwica as discussed in Section 1.5 Jand portrayed wn Fi qe lad. 
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Both the Western and Eastern Rifts north of 
about 4°S show broad (200-500 km wide) negative gravity 
anomalies. To the south of 4°S these negative anomalies 
disappear, but along the Western Rift there are localised 
positive anomalies (Wohlenberg, 1975). The gravity data 
in this region are however very sparse. North of 2°S the 
Eastern Rift has a central positive anomaly superimposed 
on the broad negative anomaly, usually interpreted as 
evidence of intrusion of mantle-derived material (Girdler 


and Sowerbutts, 1970; Searle, 1970; Baker and Wohlenberg, 
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loys Danrvacott et val. 1972 Wohlenberg, 1975). There 
is no central positive anomaly and thus no evidence of a 
central intrusion under the Western Rift (Wohlenberg, 
1975). For no rift is there geophysical evidence of 
partial melting in the mantle south of 4°S, at which 
latitude the fracture pattern in the eastern rift 
changes from grid faulting in the north to block faulting 
in the south (Fairhead and Girdler, 1972). Girdler (1975) 
drew attention to a large negative anomaly on the Bouguer 
anomaly map of Africa (Slettene et al., 1973) as shown 
ImesiMmpliftied form in Figs 6.6a. Grirdler pointed out 
that north of 10°S the anomaly can be associated with 
rifting, but ne~could not offer an explanation for, the 
westward extension. Fiqure 6.6b shows a simplified relief 
mao of Africa and it is clear that the negative anomaly 
pattern correlates with the uplifted regions of the 
continent. This means no more than that the uplifted 
topography is tsostatical ly compensated. | The east-west 
arm of this negative Bouguer anomaly lies about eight 
degrees north of the induction anomaly. 

Figure 6.7 shows the outline of the conductive 
Zone superimposed on a simplified compilation of the 
Bouguer anomaly maps of South-West Africa (Geological 
Man of South-West Africa (showing mineral occurrences 


and aravity contours), 1963; Kleywegt, 1967), Botswana 
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Figure 6.6a Simplified Bouguer anomaly map of Africa 
(Girdler, 1975) 


Figure 6.6b eee relief map of Africa (Bailey, 
Hi: 7225) 
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Figure 6./7 


Simplified Bouguer anomaly map of South-West 
Africa, Botswana and Zambia (Geological Map of 
South-West Africa, 1963; Bouguer Anomaly Map of 
Botswana, 1974; Provisional Gravity Map of the 
Republic of Zambia, 1974). 
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(Bouguer Anomaly Map of Botswana, 1974; Reeves and 
Hutchins, 1975) and Zambia (Provisional Gravity Map of 

the Republic of Zambia, 1974). The conductive structure 
and the positive gravity anomalies in the area Shows a 
remarkable correlation. Most of Botswana shows a positive 
anomaly pattern compared with Zambia and South-West Africa 
with the exception of the coastal strip. Since Tertiary 
times most of Botswana has been going through a phase 

of depression with respect to the neighbouring areas which 
has gone through several phases of uplift (Haughton, 1969). 
The changes in relief are thus in accordance with the 
qeneral gravity pattern. Where the solid rock is exposed 
in Botswana, as for instance to the south of the south- 
eastern limit of the conductor, some of the gravity highs 
are produced by Archean greenstone belts. Reeves and 
Hutchins (1975) pointed out that a change in the “texture” 
of the anomalies at a SW-NE striking line coinciding in 
the east with the south-eastern limit of the conductive 
zone divides the gravity map of Botswana into two distinct 
arease the authors speculated that™tnis cnange may be 
associated with an abrupt change in crustal thickness or 
density or that it may represent the edge of the foreland 
to the Damaran orogenic belt. Anomalous density distribu- 
tion resulting from folding and intrusion during the 


Damara tectono-thermal event could account for some of the 


ny 

aga caveo% SER! eaampEgA 7 > aot _thmean oath 
Jeo?) Une ede Seer ames Dne eur seater | 
syiP A zr TSaHube? ah (Sti ,stdned Foust idvosh oft * ; 
Ypore Ghd of SSP Uy oe evtsizoq: ond. bas aay 7 
> +o gen”  aneeéfsa-g9 ofdasteamoy va 
Tica) Aibw Gs eqeem Meme eg < lamene 
fas (at yaa nih 90 WOR gSaKe end fata 


ot sed spdvetod, tO Mzom eaaly | 


' 
vane [Le aaoedetes att erases By he noreraiesne ts 
ntigplit't ob 95 bbbete Tbe e RbOThy sa0n) Ran ee 

| | 7 


| \ iL U6? pve 34009 nt espnets Sa) oe 
aon” oi sieges iver Te vstep 


s41 fe ivtvot ot? ef Sipnsea? WOE Ge /hawesnaine 
reer re a 
a ae yl SN BAAG nawiak Aa boouborg 396) 
eas t piconet é i.e dun Deteted covet) eritneteh 
(i0a" af 4% GM-W2 2 28 corfamoas! hal Ht rae 
tee aft Se B1o8ly fr eeeg gett $oe ers 


247 de? cid) satgedad Wien ede ) 
#4, S0R882 haa tues val eluoeae 210 

“9 Deon, bal feiomee wt ie ie tourk 

na terde $46) 0) gabe any | 

-ueiivae Vb viii 2a tue eats 

: igus err neb oil capil 6a 


ont %4 ead: et Ages 


a 
ohne 


5 aah ne 


234 


gravity anomalies to the north of this boundary. Kleywegt 
(1967) found that the north-easterly trend in the gravity 
anomalies in South-West Africa is indeed related to 
folding in the Damara eugeosyncline. In South-West 
Africa the gravity trend associated with the conductor 
cuts across the NE-SW trend of the anomalies related to 
the Damaran orogeny. Figure 6.8 shows a simplified 
isostatic anomaly map for South-West Africa as computed 
by Kleywegt (1967) for an Airy-Heiskanen model with 
a sea-level crustal thickness of 30 km. The conductive 
zone coincides with the largest isostatic anomaly in the 
area. This anomaly reaches a maximum of + 660 g.u. For 
most of its lenath in South-West Africa the conductor 
coincides with an area undergoing subsidence related to 
the subsidence in the Okavango Delta and Kalahari basin. 
The reason for this correlation between an 
induction anomaly and positive gravity anomalies is not 
Clear. dhe positive gravity anomalies could be related 
tovanomalous densities and/or crustal thinning. jn 
Chapter 4 it was pointed out that the amplitudes of 
magnetic variation fields were found to be smaller in 
Botswana than further west. This reduction in amplitude 
could be related to a thinner crust or lithosphere in 
northern Botswana, and some or all of the positive gravity 


anomalies mav be due to thinning. lit Ws worth noting that 
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Simplified isostatic anomaly map for South-West 
Africa as computed for an Airy-Heiskanen model 
with a sea-level crustal thickness of 30 km 


(Kleywegt, 1967). 


figure 6,3 
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the Limpopo Belt between the Kaapvaal and Rhodesian 
cratons is associated with a positive Bouguer anomaly 
(Smit, Hales and Gough, 1962). The coincident gravity 
and induction anomalies may thus mark a similar weak 
zone between two old stable regions. In the eastern 
part they coincide with the northern boundary of the 


Kaldharw craton of Clifford (1967). 


6.2.3 Geomagnetic depth sounding and heat flow 


Geomagnetic depth sounding across the Eastern 
Rift at the Equator (Banks and Ottev, 1974) revealed 
anomalies accountable in terms of a conductive region 
20 km below the rift floor and another 50 km deep 1090 km 
Gast Of the Rett. Soundings soutn of S°S might qive 
different results, in view of the seismological evidence 
of a change in asthenosphere properties near 2°S. 

Recently eight new heat flow values ranging 
from 55-76 mWm~° were reported from Zambia (Chapman and 
Pollack, 1975). The authors compared these values 
(shown in Fig. 6.1) with the mean for Precambrian provinces 
elsewhere (Polyak and Smirnov, 1968) and found the Zambian 
values high by about 50%. This led them to conclude that 
these results sunport the hypothesis of an incipient arm 
of the East African rift svstem (a zone of crustal 
thinning) runnina southward from Lake Tanganyika to a 


point just west of Lake Kariba as proposed by Fairhead and 
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Girdler (1969). The writer believes this comparison is 
not fully justified. The heat flow results for South 
Africa as summarized by Carte and Van Pooyen (1969) give 


means, foy the 2,600 “yr old Kaapvaal shield of 50 mWm~2 


and the 1,000 Myr old Namaqualand-Natal Belt of 60 mim=2, 
The seven Zambian sites affected by the Damaran-Katangan 
tectono-thermal event between 500 and 750 Myr ago give 
a mean of 67 mim” ©, which is not abnormally high compared 


to that of the Kaapvaal shield and Namaqualand-Natal 
Mobile Belt. 


6.2.4 Origin of the conductive zone 


Regions of high conductivity which produce 
geomagnetic induction anomalies may be related to high 
temperatures and partial melting in the upper mantle, 
as in the western U.S.A. (Gough, 19/74), or to compositional 
Conductors such as graphite Or Saline Water which are orten 
concentrated in fracture zones (Gough and Camfield, 1972; 
Garland, 1975). If the conductive zone under discussion 
is associated with the African rift system its conductivity 
could be associated with anomalously high temperature in 
the upper mantle, or with a fracture zone in the 1litho- 
sphere, or with both. No definitive CeGisiOn ais DOSSip ben 
but a fracture zone seems the likeliest explanation. From 


Kariba to the Okavango the conductor can reasonably be 
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associated with the known fractures, whose activity 
resembles that of the Western, Lake Tanigamyaka Rint. 

The absence of known fractures along the westward arm 

of the conductor does not preclude their existence 

because of the Kalahari cover, and seismic aCe Wey mates 
present. The westward arm is in line with the eastern- 
Most pare Of the Walvis Ridge (Fiidyven4 )sewhichemay. on 
recent geonhysical evidence, have a fracture-zone origin, 
Sicier ta pure trans form Taullt or some tvne tof tectonic 
process following an older line of weakness (Goslin et al., 
1974; preliminary renort from Scientists aboard Glomar 
Challenger for Leg 40 of DSDP, 1975). The conductivity 
anomaly could lie on the original line of weakness 

within the continent, along the boundary in pre-Damara 
basement rocks between the eugeosyncline and the 
miogeosyncline. It is a well-established fact that 
pre-existing cratons and orogenic belts in southern Africa 
have often controlled later tectonics (Cox, 1970; McConnell, 
No7Z-eBaiveys |9724 byte and Leonardos, 1973, Thorpe 

and Smith, 1974). The seismicity may thus trace contem- 
porary failure, under a general WNW-ESE tensile stress 
filelidy on pre-existing fractures. In addition the gravity 
data indicate the crust may be tninner along the conductor. 
The first down-faulting in the middle Zambezi probably 


occurred during early Mesozoic times (Cox, 1970; Gough and 
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Gough, 1970) along fractures in the Upper Proterozoic- 
Lower Palaeozoic mobile zone related to the tectono-thermal 
Katangan and Damaran Episodes (Clifford, 1967). The 
conductor could mark this zone of weakness in the 
lithosphere (Watterson, 1975) along which intermittent 
movement has occurred since the Upper Proterozoic. 

The possibility of a thinner crust associated 
with the conductor raises the possibility that some of the 
anomalous conductivity could be due to heat. Another 
point that has to be borne in mind is that Cainozoic 
volcanism in Africa has been limited almost exclusively 
to areas affected by the Damaran-Katangan tectono-thermal 
event (Thorpe and Smith, 1974). This affinity of volcanism 
for intercratonic areas is consistent with the higher 
hear reuUxson such areas wn relation, tor the cracons 
(Fig. 6.1) and may reflect variations in other properties 
of the lithospheric mantle such as variation in chemical 
and/or mineralogic components which favour formation and 
segregation of basaltic magma (Kesicon 4973) 

| The absence of evidence for opening of rifts in 
southern Africa can be related to global tectonics in two 
alternative ways. In terms of classical plate tectonics 
the pole of opening for the East African rift system is 
in southern Africa (Fairhead and Girdler, 1971; Baker et 


al., 1972) where opening will therefore be small. On the 
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membrane stress hypothesis of Oxburgh and Turcotte (1974) 
the rifts are caused by northward movement of the continent 
over the equatorial bulge. Again propagating fractures 


with small displacements would be found in southern Africa. 


6.3 Concluding remarks 


The two magnetometer array studies discussed in 
this thesis served as reconnaissance surveys to determine 
the suitability of this method for resistivity surveys 
of the old shields and mobile ihe in southern Africa. 
The results were different than expected, and perhaps 
raised more questions than they answered. From the 
discussion in this chapter and the previous one it is 
clear that these surveys set the stage for several 
experiments, in diverse geophysical disciplines, that 
could tead to a much better understanding of the deeper 
structural differences between the various tectonic units 


in the area aS well as the resistivity structure within 


each unit. 
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